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Abstract 



We review the mechanism of electroweak baryogenesis. The main 
focus of the review lies on the development of quantum transport 
equations from first principles in the Kadanoff-Baym framework. We 
emphasize the importance of the semi-classical force that leads to re- 
liable predictions in most cases. Besides, we discuss the status of 
electroweak baryogenesis in the light of recent electric dipole moment 
probes and collider experiments in a variety of models. 
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1 Introduction 

The goal of any baryogenesis mechanism is to explain the observed asymme- 
try between matter and anti-matter 
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1] = ~ 10 . (1) 
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In order to produce such an asymmetry dynamically, several symmetries have 
to be broken what is summarized by the Sakharov conditions pQ: Baryon 
number (B) must not be conserved. Charge conjugation (C) and charge con- 
jugation in combination with parity conjugation (CP) must not be a symme- 
try. Time reversal must not be symmetry what in the early Universe implies 
a non-equilibrium state of the plasma. Owing to the Sakharov conditions, 
baryogenesis is only possible in extensions of the Standard Model (SM). In 
particular, new sources of CP violation and sizable deviation from thermal 
equilibrium are essential for a viable baryogenesis mechanism. 

The special appeal of electroweak baryogenesis [2] (EWBG) is hereby 
that only physics of electroweak scales is involved. This makes the scenario 
in principle testable. The basic picture of electroweak baryogenesis is as 
follows: At temperatures above the electroweak scale, the electroweak gauge 
symmetry is unbroken and the Universe is filled with a hot plasma of particles 
with no net baryon number. The Universe expands and cools and eventu- 
ally the electroweak gauge symmetry is spontaneously broken via the Higgs 
mechanism. Electroweak baryogenesis can be realized if this change of phase 
proceeds by a first-order phase transition. In this case, bubbles nucleate that 
contain a plasma of broken electroweak symmetry and subsequently expand 
in the surrounding plasma with unbroken symmetry. Individual particles 
in the plasma experience the passing bubble interface because of their cou- 
plings to the Higgs field. This leads to the reflection of particles and drives 
the plasma out of equilibrium. Eventually this reflection process entails CP 
violation and an asymmetry between particles and anti-particles accumulates 
over time in front of the expanding bubble walls. Since baryon number is 
conserved up to this point, the opposite CP asymmetry accumulates inside 
the bubbles of broken plasma. Finally, baryon number is violated due to the 
sphaleron process that is only active in the unbroken phase. The sphaleron 
also provides the C violation since it couples only to left-handed particles. 
This mechanism is most efficient when the particle asymmetries diffuse deep 
into the unbroken phase where the sphaleron rate is unsuppressed [3]. The 
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Figure 1: Sketch of the electroweak baryogenesis mechanism: The Higgs bubble 
walls separate the symmetric from the broken phase. If the reflection of left-handed 
electroweak particles entails CP violation, the sphaleron process (that only is active 
in the symmetric phase) generates a net baryon number. 

mechanism is sketched in Fig. [TJ 

Under all models that provide the necessary ingredients for electroweak 
baryogenesis the minimal supersymmetric standard model (MSSM) has a 
prominent role. This is mostly due to the fact that the MSSM overcomes 
(or alleviates) many shortcoming of the SM in some regions of the parameter 
space: The hierarchy problem of the Standard Model (SM), unification of 
gauge couplings, the anomaly in the gyromagnetic moment of the muon, 
viable dark matter candidates and so on. For these reasons, the MSSM is 
also the most studied framework for electroweak baryogenesis. 

Unfortunately, electroweak baryogenesis is not so easily realized in the 
MSSM. The main reasons are that the Higgs sector is rather constrained and 
that CP violation arises in a special form. Even though a strong enough phase 
transition is possible in a small region of the parameter space (the so-called 
light stop scenario), the observed baryon asymmetry can only be explained by 
nearly mass degenerate charginos and/or neutralinos. Therefore, a reliable 
analysis of the produced baryon asymmetry has to account for flavor effects 
as e.g. flavor oscillations, resonant enhancements and transport phenomena 
that are specific to the multi-flavor case. A large part of the literature deals 
with these complications that are responsible for the large discrepancies in the 
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baryogenesis analysis between different approaches. Recently, the available 
parameter space for viable MSSM models shrunk significantly with the LHC 
results, in particular the Higgs searches. All in all, electroweak baryogenesis 
in the MSSM is technically not ruled out yet, but only possible under rather 
contrived assumptions and at the cost of additional cancellations and tunings 
(a more detailed analysis will be given in section |4T4"|) . 

The main purpose of the present review is to turn the spotlight on elec- 
troweak baryogenesis in models other than the MSSM. The emphasis is 
hereby on the following aspects: 

• Using the Schwinger-Keldysh formalism, quantum transport equations 
have been derived in the recent years from first principles in the context 
of electroweak baryogenesis. Especially, when the CP violation oper- 
ative in baryogenesis results from the semi-classical force and is not 
based on flavor mixing, all applied approximations are well justified 
and allow for robust quantitative predictions. 

• Recent LHC results marked the discovery of a Higgs-like particle with 
mass rrih — 125 GeV. If this particle is identified with the Higgs particle, 
this is most relevant for electroweak baryogenesis. The strength of the 
electroweak phase transition is tightly linked to the Higgs mass. Larger 
Higgs masses tend to weaken the phase transition and suppress the 
produced baryon asymmetry. In all models we assume a Higgs mass of 
above value in this manuscript. 

• The main motivation for new physics at the electroweak scale (and 
supersymmetry in particular) comes from the hierarchy problem. The 
discovery of the Higgs highlights this fact and rules out Higgs-less mod- 
els as e.g. Technicolor. In the last years much progress was made con- 
cerning alternative solutions to the hierarchy problem as for example 
composite Higgs models. These models typically allow for electroweak 
baryogenesis without much tuning in the Higgs sector. 

The plan of the review is as follows: In section [2] semi-classical transport 
equations are derived from first principles in the Schwinger-Keldysh formal- 
ism. The main result of this section is the Boltzmann equation (148 p that 
includes a CP-violating semi-classical force at order h. Subsequently, in sec- 
tion [3] we present how to transit from Boltzmann type transport equations 
to diffusion equations using the flow ansatz and the complete analysis of the 
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produced baryon asymmetry is illustrated. Finally, in section H] the analysis 
of the baryon asymmetry and its correlation with collider phenomenology is 
discussed in specific models. The appendix contains the remaining ingredi- 
ents of the baryogenesis calculation. This includes the characteristics of the 
phase transition and the sphaleron rate. 
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2 Transport in the Kadanoff-Baym framework 

In this section, we discuss how quantum transport equations can be de- 
rived from first principles in a QFT setting and its application to elec- 
troweak baryogenesis. Main aim of this section is to sketch and motivate 
the Schwinger-Keldysh formalism (also known as closed time path formalism 
or in-in formalism) rather than to discuss it in complete depth. The dis- 
cussion closely follows the derivation in Ref. [I] and the technical review in 
refs. [5j |6]. More details can be found in Refs. [7J EJ |9]; thermal field theory 
is covered in the books [101 E] • 

2.1 The Schwinger-Keldysh formalism 

Starting point of the Schwinger-Keldysh formalism is the observation that 
not only scattering amplitudes allow for a representation in terms of path 
integrals but also the time evolution of expectation values of operators [121 
IT3] . Consider a quantum mechanical system with coordinate q, a basis n and 
some operator O that at initial time to leads to the matrix elements 

O mn (to) = (m\6\n) . (2) 

Matrix elements evaluated at later time can be related to O mn {to) via 

O ab (ti) = (a\e lA{tl - to) \m) O mn {t Q ) (n| e -**ft-*0|&) , (3) 

n,m 

Hence, unlike scattering amplitudes the time-evolution of a matrix element 
involves the evolution of states back and forth in time. 

In the path integral formulation, the evolution of the basis states can be 
expressed as 

(nle-^^-^lb^ = J Vq e if % dtc[qA) , (4) 

with the Lagrangian C and appropriate boundary conditions. The time- 
evolution of an operator can then be represented as 

CUti)= / VqO(t 1 )e i Sv^(^) j (5 ) 

using a closed time path V that goes from t to late times and back, see 
Fig. [2J It is important to remember that the two branches of integration 
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Figure 2: The closed time path contour for a general out-of-equilibrium system 
(top) and a system in equilibrium at finite temperature (bottom). 

are independent such that in the Hamiltonian picture, operators are path 
ordered and not time ordered. 

In QFT the same route can be followed leading to path integrals along 
the closed time path. As in the quantum mechanics example above the 
evaluation of operators then leads to path ordered expectation values. This 
in turn leads to the fact that the Dyson series of time-dependent perturbation 
theory does not only involve the time-ordered Green function but also the 
anti-time ordered and unordered ones. This can be expressed efficiently by 
giving the two-point functions an additional 2x2 structure, e.g. in the case 
of a scalar field one defines 

A ++ (u,v) = A\u,v) = -t(tt\T[<f){u)<j)\v)]\tt) , 

A + -(U,V) EE A < ( U ,v) EE -i(n\ft(v)(j)(u)\n) , 
A-+(U,V) EE A>( U ,V) EE -i(n\(j)(u)tf(v)\n) , 

A (u,v) ee A\u,v) ee -i(n\T [(/){u)(f>\v)]\n) , (6) 

where T and T denotes time and anti-time ordering, respectively. Obviously 
only two of the functions are independent and the matrix A in this ± notation 
is anti-Hermitian in the sense that 



A\u,v) = -A(v,u). 



(7) 
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In many cases it is advantageous to express the two-point functions in terms 
of the spectral function A = i(A > — A < )/2 and the symmetric propagator 
T = (A > + A < )/2. For canonically normalized fields, the spectral function 
fulfills the relation 

2d Uo A(u,v)\ uo ^ vo = 5(u-v), (8) 

that follows immediately from the equal time commutation relations of the 
field 0. 

Ultimately, the matrix elements (j2J) can be used to determine the prop- 
erties of a statistical system using the density matrix p 

Tr (pd) = PmnO nm . (9) 

If the density matrix is known at initial time, operators can be evaluated at 
later times using the path integral representation of O mn as outlined above. 
In principle all information about the system can then be inferred from the 
density matrix. 

An alternative way of proceeding is to consider a closed system of n- 
point functions and to impose the initial conditions on the n-point functions 
rather than the density matrix. In complete analogy to QFT calculations, the 
Schwinger-Dyson equations can be derived from the 2PI effective action [TJ] 
in the non-equilibrium setup. Formally, the equation is the same, namely 

/ d 4 w(D + m 2 +U(u,w))A (w, v) = 8{u — v) , (10) 

where IT denotes the self-energy. In a specific model, the self-energy II can 
be expressed perturbatively in terms of the interactions and the two-point 
functions of the system. This allows to determine the two-point functions 
at all times consistently without resorting to initial conditions in terms of a 
density matrix 0. 

However, even though these equations are formally the same as the Schwinger- 
Dyson equations, the two-point functions are understood to have the addi- 
tional 2x2 structure mentioned before. Besides, in many cases statistical 

1 Strictly speaking, the Schwinger Dyson equation in the 2PI formalism allows only for 
Gaussian initial conditions. More general initial conditions require the use of the nPI 
formalism or similar techniques JTSJ El [TBI H3 HH] • I n the present context, this problem is 
not of relevance, since one applies the limit to — > — oo and hence thermal initial conditions. 
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systems are not isotropic or homogeneous such that the two-point functions 
A and the self-energy II do not only depend on the relative coordinate (u — v) 
but explicitly on both coordinates u and v separately. This feature is par- 
ticularly bothersome if the two-point functions are transformed into Fourier 
space. Usually Feynman calculus is particularly simple in Fourier space since 
the convolutions in coordinate space turn into conventional products 

J dyA(x - y) B(y - z) ^ A(p) ■ B(p) . (11) 

However, if a dependence on the average coordinate remains, convolutions 
turn into Moyal star products 

J dyA(x,y)B(y,z)^A(p,X)*B(p,X). (12) 

Here A(p, X) denotes the Fourier transform with respect to the relative co- 
ordinate r = (x — y) for fixed central coordinate X = (x + y)/2 

A(p, X) = J d 4 r A(X + r/2, X - r/2) e ir ' p , (13) 

and the Moyal star product is defined using the diamond operator 

o = I (K^* - *dx1f P ) , (14) 

by 

A(p, X) * B(p, X) = A(p, X) e-^B(p, X) . (15) 

This representation of two-point functions is called Wigner space and allows 
for an interpretation in terms of a semi-classical phase space. One particu- 
larly simple application of this formalism is QFT at finite temperature what 
we discuss next. 

QFT at finite temperature 

The density matrix at finite temperature is given by the Hamiltonian H and 
the temperature T = as 

p = cxp(-H/3) . (16) 
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The partition function of this system 

Z = Trp= I dq (q\e- 6 %) , (17) 



can be represented by extending the closed time path into the imaginary 
time direction (see Fig. [2]) and imposing periodic (anti-periodic) boundary 
conditions for bosonic (fermionic) fields. For the two-point functions, the 
periodic boundary conditions turn into the Kubo-Martin-Schwinger relation 



A> (u,v)\ uo _ VQ=t = A < (u,v) luo _ vo=t+i ^ 

F ' T 4 A>(k) = eMk P)A < (k). (18) 



In combination with the spectral sum rule (jSJ) 

J^2poA(p) = l, (19) 

this yields in equilibrium for a free field 

A(p) = 7i5(p 2 - m 2 ) sign(po) , 

Tip) = -7ri5(p 2 -m 2 )[2n(\p \) + l] , (20) 

or equivalent ly 

A < = —ixi 5(p 2 — m 2 ) sign(po) n(po) , 

A > = — ixi 5{p 2 — m 2 ) sign(p ) (n(p ) + 1) . (21) 

Here we recover the Bose-Einstein particle distribution function 

= exp( ^) _ i ■ < 22 > 

For particle species that are weakly interacting and close to equilibrium, 
the spectral function A is approximately still given by a 8 function and the 
corresponding component of the plasma can be described by quasi-particles. 
The particle distribution function n(X,p) is then encoded in the symmetric 
propagator J 7 or the Wightman functions A < ' > . 
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Kadanoff-Baym equations 

It is not surprising that the Wightman functions A < and A > encode the 
particle densities in the plasma. After all they represent the particle number 
operators. This indicates a way to derive quantum transport equations from 
first principles: The Schwinger- Dyson equations ( jTUi) in Wigner space (that 
are also called Kadanoff-Baym equations [19|) 

(pa _ m 2 + n(p, X)) * A(p, X) = 1 , (23) 

have to be solved with appropriate boundary conditions. In components this 
equation can be brought to the form [51 [6] 

(p 2 - m 2 - U h ) * A<'> - n<'> * A h = coll. , (24) 
where we introduced the collision term 

coll. = l - (IP * A< - n< * A>) , (25) 
the Hermitian part of the Green function 

A^ = A'-^(A < + A > ), (26) 

and analogous definitions for the self-energy II. Once the Wightman func- 
tions are known, the particle distribution functions can be read off at late 
times when the system is again close to equilibrium. According to ( 12T|) one 
finds 

n(X^p) = 4i / ^A<, (27) 
Jp >0 111 

l + n(X»,p) = Ai ! ^A<. (28) 
Jpo<0 271 



Using appropriate boundary conditions, the equations ( l24j) can be readily 
applied to the problem of electroweak baryogenesis. Initially the system is 
close to equilibrium and during baryogenesis driven out of equilibrium. In 
the case of electroweak baryogenesis this stems from the bubbles of Higgs 
vacuum expectation value that give rise to a space-time dependent mass 
term m(X). The terms on the left-hand side describe the forces that act 
on the particles and also the diffusion of the particle densities away from 



2 TRANSPORT IN THE KADANOFF-BAYM FRAMEWORK 13 

the wall. The term on the right-hand side (that is called collision term) 
represents the interactions that drive the system to kinematic and chemical 
equilibrium. The particle densities of the species under consideration can 
then be read off from the Wightman functions at late times after the phase 
transition is completed. 

2.2 Approximation schemes 

In order to make the system of equations more manageable, several 
approximations can be applied that we discuss in this subsection. In the 
context of electroweak baryogenesis the following approximations are usually 
employed 0: 

• Gradient expansion: 

If the background depends only weakly on space and time coordinates, 
an expansion of the Moyal star products in the diamond operator can be 
performed. Naively, this is a good expansion for electroweak baryogen- 
esis since the background is only slowly varying in units of the typical 
momentum scale. To be specific, in the MSSM the thickness of the 
Higgs bubble wall is typically of order £ w ~ 20 — 30 T -1 . At the same 
time, a typical particle in the plasma has a momentum of order p ~ T. 
Hence, the diamond operator comes with a factor o ~ (£ W T) _1 <C 1. 

• Fluid approximation: 

The plasma is assumed to be close to equilibrium. In particular, it is 
assumed that two-to-two scatterings (or other interactions that do not 
change particle numbers) are fast such that the plasma is well described 
by the local velocity of the different components of the plasma, the local 
temperatures and the chemical potentials. The particle distribution 
functions can then be parametrized as 

1 

n — 

e (n Al pM+ A1 )/3 _j_ j ' 

where u^, (3 and fi are space-time dependent and denote the four- 
velocity, the inverse temperature and the chemical potential of the 

2 See Ref. |2Q| for a similar discussion. 



(29) 
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components of the plasma. By taking different moments of the trans- 
port equations, the equation of motion for these quantities can then be 
derived (this is exemplified in section [3]). 

• Weak coupling: 

Far away from the source of non-equilibrium the system will reach 
its chemical equilibrium via interactions that change particle numbers. 
These interactions are assumed to be slow such that an expansion in 
the according coupling constants can be performed. 

In light of these assumptions, one can then simplify the Kadanoff-Baym equa- 
tions. As a word of caution, notice that the validity of these approximations 
is not always guaranteed. The prime example are flavor oscillations where the 
fluid approximation can fail [20]. We will comment on this issue in sec. 12.51 
We will see shortly that the main source that drives the system out-of- 
equilibrium and induces CP violation arises from a kinematic effect that even 
persist in the limit of vanishing interactions. The deviations from equilibrium 
are then suppressed by e w ~ (l m T) _1 while the self-energy is suppressed by 
coupling constants and loop factors, e co u ~ g 2 /4ir. In particular, the collision 
term vanishes in equilibrium but also has an explicit factor e co u from the 
self-energy. Hence one can neglect the higher gradients in the Moyal star 
product of the collision term and write 

coil. ~ n> (p, x) a< (p, x) - n< (p, x) a> (p, x) . (30) 

Furthermore, the terms involving the self-energy on the left-hand side of 
the Kadanoff-Baym equation ff24|) mostly affect the shape of the spectral 
function. The term involving U h renormalizes the mass term while the term 
involving Il < ' > leads to a broadening of the spectral function [51 [6]. These 
terms will also be neglected in the following such that the Kadanoff-Baym 
equations read 

(p 2 -m 2 ) * A<'> = n > (p,X)A < (p,X)-n < (p,X)A > (p,X). (31) 
2.3 One bosonic flavor 

For a system with only one bosonic degree of freedom, the Wightman func- 
tions are purely imaginary and one can immediately split the Kadanoff-Baym 
equations into a real 

(p 2 -m 2 )cos(o)A < - > (p,X) = 0, (32) 
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and imaginary part 

(p 2 - m 2 ) sin(o) A < - > (p, X) = coll. . (33) 

The real part determines the spectral function and is usually called constraint 
equation while the imaginary part describes the variation of the particle dis- 
tribution functions due to the background and is called kinetic equation. 
To order 0(o 2 ) these equations are solved by the ansatz 

A< = 2n5(p 2 - m 2 )sign(p ) n(p, X) , (34) 

where the particle distribution function now fulfills the equation 

2n5(p 2 - m 2 ) (2p^d fl n(p, X) + d fl m 2 (X)d p ^n(p, X)) = coll. (35) 

This equations allow for a simple semi-classical interpretation: Imagine a 
particle with a space-dependent mass m 2 (z) arising from the Higgs bubble 
and a fixed four-momentum p M in front on the bubble wall. If the particle 
passes the wall, its mass changes. If the semi-classical particle is on-shell 
on both sides of the wall, it has to change its four-momentum and the sym- 
metries of the problem dictate that this change arises in p z . This reasoning 
leads to the relation p 2 in + m 2 n = p 2 out + rr? out and the approaching particle 
perceives the change in mass similar to a potential barrier. In particular, very 
soft particles cannot fulfill the on-shell condition inside the bubble and are 
reflected by the bubble wall. If this picture is generalized to a distribution 
of particles, n(p, X), and a smoothly changing mass profile m(X), this leads 
to the statement 

p^n^X) = -d,m 2 (X)d p n(p,X) , (36) 

which is eq. (135]) in the absence of interactions. In the language of Boltzmann 
equations, the change in mass leads to a kinematic effect that exerts a force 
on the particles in the plasma. This effect is purely classical in the sense that 
it will not be suppressed in the limit h — > 0. For electroweak baryogenesis, 
this effect is interesting since, as we will see, in case of fermions and/or several 
flavors the kinematic forces can entail CP violation (to first order in h). 

Before we do so, let us comment on some additional features of eq. fl35l) 
and its solution. First, notice that if the wall is at rest relative to the plasma, 
the force is absent. In the wall frame the mass depends only on the spatial 
coordinates, m(z), while in the plasma frame the equilibrium distribution 
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function depends only on the energy, n(po). If these two frames coincide the 
force term d^m 2 (X)d p n(p) vanishes and the equilibrium solution (with the 
space-time dependent mass) solves fl35|) everywhere. In terms of particles, 
soft particles are reflected, n{p) = n(—p), while hard particles replace hard 
particles on the other side of the bubble wall. During electroweak baryoge- 
nesis, deviations from equilibrium are hence additionally suppressed by the 
wall velocity v w in case it is substantially smaller than the speed of light. 

Next, notice that the effect persists even in the limit of vanishing in- 
teractions. Once the wall is moving, the soft particles are still reflected, 
n{p) = n(—p), but this is not consistent with the boundary conditions of a 
plasma moving towards the bubble wall. Also behind the wall the plasma is 
not in equilibrium. So interactions are essential to establish equilibrium far 
from the wall but are not so important to generate the out-of-equilibrium 
situation in the present context. 

Finally, notice that as long as the effect from the wall can be expressed 
as a force 

p^n(p,X) 5(p 2 - m 2 ) = m(X)F»{X)d Pfl n(p,X) 5(p 2 - m 2 ) , (37) 
the four-current 

J» = J d^pp^niX.p) 5{p 2 - m 2 ) , (38) 

is conserved 

«V M = 0. (39) 

This supports the picture that the effect is kinematic and neither are particles 
created nor destroyed in the process. Of course, including particle number 
changing interactions from the collision term modifies this conservation law. 
On the other hand, energy-momentum 



= I (Fpp»p u n(X, p) 5(p 2 - m 2 ) , (40) 

is not conserved 

<9 M T^ = / <fp m F v n(X, p) S(p 2 - m 2 ) ^ , (1 1 ) 



due to the latent heat that is released during the phase transition from the 
Higgs sector into the plasma. But interactions preserve the (total) energy- 
momentum tensor 

|^ coll. = (42) 
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and do not modify this relation (when summed over all species). Ideally, 
any approximation to the transport equations that is applied subsequently 
should respect these laws. 

2.4 One fermionic flavor 

In case of a system with one fermionic flavor, the derivation of the Kadanoff- 
Baym equation parallels the bosonic case. The equation in correspondence 
to (IHT1) yields in this case 

ty-P L m(X)-P R m*(X))*S < (p 1 X) = coll., (43) 

where S K denotes the fermionic Wightman function. All subleading terms 
are already neglected and we introduced a complex, space-time dependent 
mass. Unlike the bosonic case, this equation cannot be simply split into 
constraint and kinetic equation because the Dirac operator as well as the 
Green function S < contain a spinor structure. In the following we assume 
that the change of the mass is aligned with the momentum of the particle 
(both in ^-direction in the following) what makes the problem effectively 1 + 1 
dimensional. In case these two directions are not aligned, this situation can 
be achieved by a suitable Lorentz boost [2Tj . 

The spinor structure can then be partially decoupled by observing that 
the Dirac operator commutes with the following spin operator 

S z = 7o7375 oc 7i7 2 . (44) 

Using the projectors P s = ~(1 + sS z ), the Dirac operator can be brought to 
block diagonal form. The block that encodes the vector and axial currents 
can then be parametrized as 

S< = J2PsS<, S< = P s [ l0 g s + l3 gl + g s 1+l5 g s 2 }. (45) 

s=± 

In this notation the s-even (odd) parts of go encode the vector density (ax- 
ial z-current), #3 encode the vector ^-current (axial density) and gi/ 2 the 
scalar/pseudo-scalar (z-spin densities). 

In the gradient expansion the spinor structure of the Kadanoff-Baym 
equations can be decoupled [I] what leads to the following constraint and 
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kinetic equations for g 

(k s d z -±\m 2 \%.-^{\m\*eyd k ^g' Q = coll., (46) 

where the mass term was parametrized as m(z) = \m(z)\exp(9(z)). So the 
function g^ allows again for the ansatz 

g* oc 6(k 2 - u 2 )n s , u 2 = k 2 z + \m\ 2 + f \m\ 2 6' , (47) 

with 

(k z d z - \\m 2 \'d kz - ^{\ m \^')'d k )j n s = coll. . (48) 

The additional CP-violating force in this equation leads to CP-violating de- 
viations from equilibrium in the axial ^-current. The analogous equation for 
#3 shows no dependence on the shift in phase 8'. In total, no particles are 
produced or destroyed. Still, particles with different spins perceive different 
potential barriers and are reflected differently by the wall. The spin of the 
particles is hereby conserved while the chirality is not. 

If the wall is at rest, n does not depend on k z and away from the wall the 
particle distribution functions are in their local equilibrium form. The on- 
shell condition is still different for particles with different spins such that the 
two-point functions and also the axial current J| depends on the change of 
phase 6' in the wall. Since the solution is consistent with the KMS relation, 
including interactions does not change this picture [B]. Only if the wall 
velocity is nonzero the CP violation can diffuse into the symmetric phase 
and give rise to sizable baryogenesis. 

The equation (148)) is the central relation for electroweak baryogenesis 
with one flavor. The forces on the left-hand side of the equation encode 
how the plasma is driven out-of-equilibrium and how CP violation manifests 
itself in the particle densities. The kinetic term in combination with the 
collision terms dictate how the particle densities diffuse away from the wall. 
The collision terms also determine how the asymmetries are communicated 
to the other particle species and finally the weak sphaleron. The complete 
electroweak baryogenesis calculation in a toy model is sketched in sec. [31 
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2.5 Several flavors 

If several flavors are considered, additional complications arise. The diamond 
operator comes with a factor h and for one flavor the constraint equation is 
in leading order algebraic. Besides, the kinetic equation has an overall factor 
H and is in leading order a classical transport equation. For several flavors, 
the leading order of the kinetic equation (in the case of bosons) becomes 

2P<9 M A< + i[m 2 , A<] - ~ {m 2 ', <9 fcz A<} = coll. . (49) 

The first two terms of this equation describe flavor oscillations with a fre- 
quency oo ~ Am 2 /k z oc 1/h, while the third term gives forces similar to what 
was found in the one flavor case. The Wightman function does encode in 
the case of several flavors not only semi-classical particle distribution func- 
tions but also coherent superpositions of different mass eigenstates. Even 
though the Wightman function is diagonal in mass eigenbasis far away from 
the wall, the forces induce off-diagonal entries that participate in the flavor 
oscillations. This mechanism gives rise to new sources of CP violation. In 
particular, this effect arises already in leading order in the kinetic equation. 
In comparison, the semi-classical force found for one flavor contains one more 
gradient (and hence one more factor h). On one hand, this indicates that the 
flavor mixing effects can be enhanced relative to the semi-classical force. On 
the other hand, if the oscillation is rather fast, this suppresses the efficient 
population of any off-diagonal densities. So it is not a priori clear if the 
CP violation stemming from mixing or the one from the semi-classical force 
dominates the produced baryon asymmetry. 

For completeness, we quote the kinetic equation for fermions with several 
flavors as derived in [22] up to second order in gradients. In this case, it 
is more appropriate to parametrize the two-point functions in terms of left- 
handed and right-handed densities. The equation of the right-handed density 



2 TRANSPORT IN THE KADANOFF-BAYM FRAMEWORK 



20 



reads 



i 1 
k z d z g R + - [m^m,g R ] - - {(m ] m)\ d kz g R } 

+ (m v mg R + g R m)m') - — (m v g L m + rrJg L m') 
- ^ [(rrJm)", d 2 kz g R ] + ^— [m v m', d kz g R ] 

+ K m, " m *'(f)- Sfe (S) m,m ") 

~v d k , — \ m- rrvdu, — 



- \ in O u . | —J m - m'0 kz | in 
= coll.. (50) 

The corresponding equation for the left-handed density is obtained by the 
replacements 

9r^ 9l m^m ] . (51) 

Notice that this equation does not explicitly depend on the spin quantum 
number s and we dropped the superscript. The dependence on s appears 
again when the functions are rewritten in the previous notation via 

9l = 9o-sg s 3 , g s R = g s + sg s 3 , (52) 

and the lowest order relation k z g 3 = k g . Also notice that this kinetic 
equation does not explicitly depend on the energy k . Hence, the transport 
equations for the particle distribution functions can be obtained be integra- 
tion without knowledge of the spectral function. 

The second term in (1501 induces flavor oscillations while the remaining 
term of first order in gradients are analog to the classical forces in the one 
flavor case. These terms source the off-diagonal entries (in flavor space) of the 
Wightman function and contain new sources of CP violation as in the bosonic 
system with several flavors. The last two terms reproduce the semi-classical 
force known from the one flavor case. 



Application to the MSSM 



The main application of the equation (|50|) is chargino (or neutralino) driven 
electroweak baryogenesis in the MSSM. In this framework, the semi-classical 



2 TRANSPORT IN THE KADANOFF-BAYM FRAMEWORK 



21 



force that drives electroweak baryogenesis in the one flavor scheme is insuf- 
ficient to account for the observed baryon asymmetry This is mainly due 
to a weak phase transition and rather strict constraints from EDM measure- 
ments. Hence electroweak baryogenesis in the MSSM has to rely on flavor 
mixing effects that nominally are suppressed by one less order in the gradient 
expansion. 

Unfortunately, some of the usual assumptions used in electroweak baryo- 
genesis calculations potentially break down in the case of CP violation stem- 
ming from flavor mixing as is discussed in detail in [20]. The oscillation 
frequency is of order t~J, ~ Am 2 /p where the Am 2 denotes the difference of 
the mass eigenvalues squared. For soft particles this leads to fast oscillatory 
behavior and numerically this fast oscillation suppresses the relevance of the 
off-diagonal entries. On the other hand, flavor oscillations are important for 
the new CP-violating terms that arise in the kinetic equations ( 150]) beyond 
the semi-classical force [22J. 

If the oscillations are generally assumed to be faster than the background 
gradients, r osc <C £ w , the system is in the adiabatic regime [20] [23] . In 
the case of the MSSM this seems reasonable since the bubble wall is rather 
thick, £ W T = 10 — 20, and the charginos are never mass degenerate in the 
wall. Hence, the assumption r osc <C £ w should be valid for a typical particle 
in the plasma with p ~ T. In this regime the flow ansatz (129]) (including a 
collective oscillation) seems reasonable. Besides, backreactions from the off- 
diagonal densities on the diagonal ones are smalo and can be neglected. This 
is the route followed in [21] . Unfortunately, the resulting baryon asymmetry 
is too small to be simultaneously in accord with EDM constraints and the 
observed baryon asymmetry (a more extensive account of these results is 
given in sec. 14.41) . 

A first study that does not rely on the assumption of fast oscillations 
was presented in [201 123] for a toy model. In this regime, the interplay of 
off-diagonal and diagonal parts in flavor space is more involved what can 
lead to a parametric enhancement of CP violation in the diagonal particle 
densities. In a bosonic toy model, the modes that are most affected by CP 
violation are the ones where the oscillation frequency is comparable to the 
background gradients, t osc ~ £ w . As argued before, in the MSSM these 
particles are rather hard and this leads potentially to a suppression since 

3 Nominally they are second order in gradients and compatible with the semi-classical 
force terms. 
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these hard modes are not very abundant in the plasma. Still, it might turn 
out that these modes contribute more to the CP-violating particle densities 
than the bulk of particles in the adiabatic regime. To settle this issue would 
require an analysis along the lines of [201 123] in a fermionic system (namely 
the chargino sector of the MSSM) which is a daunting task. 

2.6 Other approaches 

In this section we briefly discuss to what extent the approach presented in 
the last section is consistent with other methods found in the literature. In 
particular we discuss the semi-classical force in the WKB approximation and 
the mass insertion formalism. 

Semi-classical force in the WKB approximation 

Historically, the semi-classical force was initially found in the WKB approxi- 
mation [251 ESI EH M\ and subsequently applied to the MSSM [2H1 EDI Ell E2] • 
The derivation is a little less clean than the one in the Kadanoff-Baym frame- 
work. For example, it relies on the quasi-particle picture what is a stronger 
requirement than the mere gradient expansion used in the KB approach. 

The derivation goes as follows: Assume again one fermionic particle 
species with a space-time dependent complex mass term, m = \m\e' ie . The 
corresponding Lagrangian is 

£ = $(i@ - P L m- P R m*)ip . (53) 

Using a local axial transformation, the Lagrangian can be brought to a form 
where the mass term is real, but an axial gauge field appears 

C = $ (ip + 75 ^ - m)V> , (54) 

where = \dpQ. Solving the Dirac equation then leads to the dispersion 
relation of the quasi particles. In the wall frame one finds [261 12Z] 

E 2 =pl + ^p 2 z + m*±Z z y . (55) 

The different signs denote hereby the spin in z-direction in the frame with 
vanishing p± analogue to the construction in the Kadanoff-Baym approach 
(144|) . The group velocity of the particle is given by 

dE , . 
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and energy conservation gives the constraint 

dE dE 

£ = ° = z a7 + ^' (57) 

and hence p z = —d z E. From these relations the Boltzmann equation can be 
derived 

cLti 

— = d t n + zd z n + p z d Pz n = coll. . (58) 

Notice that the relation p z = —d z E ensures that for a static wall the equilib- 
rium particle distribution function (that in this case only depends on energy 
in the wall frame) is a solution to the Boltzmann equation. 

Let us compare this result with our findings in the Kadanoff-Baym ap- 
proach. In the 1+1 dimensional case and for small gradients one finds 

E 2 = (y/p> + m*±Z z y 

~ pl + m 2 ±2EZ Z (59) 



Comparing with (1461) we see that the force in the WKB approximation is 
smaller by a factor m 2 /E 2 what is close to unity for non-relativistic parti- 
cles. So the result is in rough agreement with the ones later obtained in the 
Kadanoff-Baym framework. However, the CP-violating term arises through 
an (axial) gauge transformation what initially lead to some discussion in the 
literature if this effect is physical. This issue can be resolved by distinguish- 
ing between canonical and physical momenta [31J. This careful analysis also 
recovers the factor m 2 /E 2 and is then in full agreement with the result from 
the Kadanoff-Baym framework. 

In conclusion, the derivation of the leading order effect in the Kadanoff- 
Baym framework agrees with the one in the WKB approximation for one 
fermionic flavor. Nevertheless, the Kadanoff-Baym framework overcame some 
shortcomings of the semi-classical analysis. First, above ambiguity involving 
the canonical and physical momenta never arises. Second, the Kadanoff- 
Baym framework does not assume quasi-particle states from the start. The 
quasi-particle properties are rather a consequence of the constraint equations 
to the lowest orders in the gradient expansion. 



Mass insertion formalism 



Another approach to CP- violating sources in transport equations is the mass 
insertion formalism [331 Ei E3 ESI EH ESI ESI SHI El S2J H31 IB] • The for- 
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malism has compared the full-fledged Kadanoff-Baym treatment the advan- 
tage that it is perturbative what makes even calculations with several flavors 
straight forward. The main application of this formalism is hence electroweak 
baryogenesis in the MSSM. 

The main idea is to treat the mass term as an interaction and expand the 
Kadanoff-Baym equations around a plasma with vanishing mass. Formally, 
the fermionic equivalent of equation (12"5|) 

(jf> - P L m - P R rrJ - E(p, X)) * S(p, X) = 1 , (60) 

is solved perturbatively (neglecting the terms £ arising from 'real' interac- 
tions) 

j, * = [P L m + P R m)) * S {0) , (61) 

where S* denotes the equilibrium solution of a massless particle. 

On general grounds this formalism gives rise to several objections |43] : 

• In the case of one flavor, the main effect comes from a shift in the 
dispersion relation. This effect can only correctly be accounted for if 
the Kadanoff-Baym equations are solved. In the perturbative picture, 
the Kadanoff-Baym equations resum an infinite set of diagrams. Even 
worse, if the operator j> * in (I6TI) is inverted one encounters divergences 
that have to be dealt with. As a simple example, consider the following 
equation that mimics the constraint equation 

(a; - a - Aa)f{x) = , (62) 

with the solution f(a) oc 8{x — a — Aa). If the equation is expanded in 
Aa, one finds p°'(a) oc 8(x — a) and f^(a) = Aa f^/(x — a), what 
is not well defined. The correct behavior can in principle be recovered 
when one identifies S(x — a)/(x — a) — > —8'{x — a). However, in the 
literature on electroweak baryogenesis the problem is usually avoided 
by introducing a finite width in the spectral function. Potentially, this 
leads to an overestimation of the effect. Without expanding in Aa, the 
result is manifestly finite. 

• By construction, the resulting Wightman function is local and hence 
does not contain any transport. To overcome this problem, the re- 
sulting deviation from equilibrium is interpreted as a source term and 
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subsequently inserted into a transport equation to make diffusion pos- 
sible. In the literature different proposals exist how this has to be done, 
the most plausible being the use of Fick's law [42] • 

• Flavor oscillations are not correctly reproduced in the studies based on 
the mass insertion formalism. 

• Once the source is inserted into the (classical) transport equations, a 
basis choice has to be made. The observation is that the CP-violating 
source vanishes in the mass eigenbasis and the interaction eigenbasis 
is used. On the other hand semi-classical quasi-particles propagate as 
mass eigenstates what makes this choice questionable. The transport 
equations obtained in the Kadanoff-Baym framework are in principle 
basis independent @. 

In ref. [42] a refined version of the mass insertion formalism was presented. 
The mass was hereby expanded around a fixed point 

m(X) = m(X°) + (X„ - Xl)d^m . (63) 

The derivative term was again treated as an interaction while the mass term 
was incorporated in the lowest order solution S^°\ This overcame some of 
the problems listed above but also reduced the predicted baryon asymme- 
try by one order of magnitude. In this partially resummed form, the main 
differences between the mass insertion formalism and the Kadanoff-Baym 
equations seem to be how transport is implemented and the neglect of fla- 
vor oscillations. While the Boltzmann type equations arise naturally in the 
Kadanoff-Baym equations, the mass insertion formalism still requires to use 
Fick's law or some other classical input to describe transport. 

A quantitative comparison between the different approaches in case of 
the MSSM is given in section 14.41 



4 However, in practice also a basis choice is often made in the Kadanoff-Baym framework 
when the particle densities are coupled to other species, see e.g. [24] • So the problem is 
for the most part only postponed. 
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3 Electroweak baryogenesis: A toy model 

In this section, we connect the analysis of CP-violating particle densities with 
the explicit calculation of the baryon asymmetry. Namely, we discuss how 
to transit from Boltzmann to diffusion equations (mostly in the case without 
flavor mixing). Finally we exemplify the complete calculation in a toy model. 
Some ingredients, as e.g. the weak sphaleron rate and the characteristics of 
the phase transition are covered in the appendices. 

3.1 From Boltzmann to diffusion equations 

To solve the partial differential equations (I48p or ( I50p is rather demanding 
without using further approximations. In the following, we discuss only the 
diffusion equations in models without flavor mixing, where the semi-classical 
force is the dominant source of CP violation. 

Consider a Boltzmann type equation in the wall frame of the form 



To simplify these partial differential equations further, often the so-called 
flow ansatz is used. The underlying assumption is that equilibration involves 
different time scales [IHIHT]. When out of equilibrium, the system establishes 
after a short time kinetic equilibrium by decoherence effects and scattering 
processes. After this phase, the particle distribution functions of individual 
species are approximately of the form 



where denotes the plasma four-velocity, T the temperature and /z the 
chemical potential. At intermediate times, these quantities are still space- 
time dependent. Only at later times, the temperature and the four- velocity 
of the different species equilibrate to each other and the chemical potentials 
approach an equilibrium consistent with the conserved charges of the system. 
Similarly, in electroweak baryogenesis the flow ansatz is fulfilled reasonably 
well everywhere while the correct equilibrium is only attained away from 
the wall. Furthermore, it is usually also a good assumption in electroweak 
baryogenesis to use the same temperature for different species. This is owed 
to the special structure of the CP-violating source^. 

5 In the calculation of the wall velocity this would be a poor approximation [48 . 



p z d z n(p) + m F z d Pz n{p) = coll. . 



(64) 




(65) 
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Before we solve the Boltzmann equations using this ansatz, we discuss 
in a little bit more detail the connection to the Kadanoff-Baym equations 
of section [2j In the Kadanoff-Baym equations, the distribution functions for 
anti-particles are given by the negative frequency part using the identification 

n(p°) = -n(-p°) ± 1 . (66) 

Hence anti-particles come in the flow ansatz (I65I) with the same four-velocity 
and temperature but opposite chemical potential, as it should be. Of course, 
in the presence of CP violation small deviations between the chemical poten- 
tials and velocities of particles and anti-particles can arise. Another impor- 
tant point is how to connect the CP-violating force to the system of Boltz- 
mann equations. The Boltzmann equations do not contain the full Dirac 
structure of the Kadanoff-Baym approach but only parametrize the system 
by four densities of (pseudo-) particles. Typically these are chosen to be 
left-/right-chiral particles/anti-particles. In contrast in the Kadanoff-Baym 
approach, spin is a conserved quantum number. In order to translate the 
semi-classical force ( HHjl into the chirality basis, the force is transformed into 

P !| 21/ , si g n (P Z ) /, 2.fl/y fR v\ 
mb z ~ ~\m | ± (\m \v ) , (67) 

where opposite signs apply for left /right-chiral densities and particles/anti- 
particles respectively. Strictly speaking this identification is only true for 
highly-relativistic particles, but we will see below that it reproduces (in lead- 
ing order in wall velocity) the correct deviation from equilibrium in terms of 
vector and axial currents. 

Using the flow ansatz, different moments of the transport equation 
can then be taken in order to reduce the Boltzmann type equation to a 
diffusion type equation. This leads to the relations 

(Pz) ^ + (pI) < + (mF z ) u z = (coll.) , 
(pI) + (pi) < + ( Pz mF z ) u z = (p z coll.) , (68) 

where u z is in leading order given by the flow of the background (that equals 
the wall velocity far away from the wall). We used the fact that the flow 
term and the force fulfill the relations du/dp z = p z /u and duj/dz = mF z /u 
what ensures that the two derivative terms acting on the energy u cancel 
each other 0. 

In principle there arises an additional term from derivatives acting on the term sign(p 2 ) 
but these turn out to be negligible [49j . 
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The moments are usually defined as 




with a normalization to a fermionic massless degree of freedom in equilibrium 

(70) 

m=/j,=u z =0 

In the following we linearize the system in the chemical potentials and the 
flow velocities. In leading order certain moments are then related by Lorentz 
boosts, e.g. (p z ) ~ — u z k where k = (u) denotes the statistical factor that is 
1 (2) for massless fermions (bosons) for a plasma at rest. Furthermore (p 2 z ) 
is in leading order 1/3 of the pressure in the plasma and (pi) ~ — 3u z (p 2 z uj). 

Next consider the collision terms. The collision integral in the second 
equation is dominated by elastic scatterings 

(p 2 coll.) ~-T ela (u-u), (71) 

(notice that T ela has dimension three according to this definition). The func- 
tion u denotes the flow velocity of the background the species mostly scatters 
with and it is often assumed that this is given by the wall velocity that de- 
scribes the flow far away from the wall, u v w . Notice that this approxima- 
tion is in principle not consistent with the arguments of energy-momentum 
conservation discussed in sec. 12.31 Still, as long as the background represents 
a large number of degrees of freedom, this approximation is reasonable. 

The collision term in the first equation encodes the particle changing 
interactions. These are of the form 

(co\l.)~T inela J2wi, (72) 

i 

with q some integer constants and the subscript i labels the species of the 
chemical potentials /j*. One of these interactions constitutes the sphaleron 
rate that finally biases the baryon number. Both sphaleron rates, strong and 
electroweak, are non-perturbative and cannot be recovered from the collision 
term as given in ( 125]) . They have to be added by hand to the network of 
transport equations. 



N 



d 3 p 



drif 
d/u, 
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Finally, consider the forces in the diffusion equation. The CP-conserving 
force drives the flow of the particles and anti-particles equally away from 
equilibrium 

(mF z )u z ~ ~|m 2 |'(lK = <V (73) 

but does not have a large impact on the chemical potentials. The CP- 
violating force on the other hand contributes mostly to the equation involving 
the chemical potential 

(p z m F z ) u z ~ ~(|m 2 |0')'(N/wK = S u . (74) 

In particular, the CP-violating force comes with different signs for the left- 
and right-chiral fields such that it has only an impact on the axial current 
as found in the Kadanoff-Baym approach. Besides, it vanishes explicitly for 
static walls. 

This system of equations can be brought to the form of a diffusion equa- 
tion by neglecting terms that are second order in the velocities in the second 
equation. This gives 

(u-u)^^((p 2 z )^ + S u ) . (75) 

Neglecting derivatives acting on the averages and using this in f|68|) yields 

Dfi" + v w Kfj! + S^ + S D = (coll.) , (76) 
where we defined the diffusion constant 

and the CP-violating source of the form Sd = S' u /T ela . However, there is no 
need for these additional approximations and the linear differential equations 
(|68|) can be easily solved numerically. 

In conclusion, the system of transport equations can after linearization 
in the velocities and the chemical potentials be brought to the form 

9 z AJ^ + ^rr /a c^ = 0, (78) 

A,/3 
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and 

d z AT* z + Yi a (u a - v w ) = S a , (79) 

where the indices a and /3 run over all particle species and chiralities and the 
index A over all interactions. The term A J z denotes the current of particles 
minus antiparticles and the expression d z AJ z represents all three terms on 
the left hand side of the first equation of (1681) . Likewise, the term AT ZZ 
denotes the zz component of the energy momentum tensor of the particles 
minus antiparticles and the expression d z AT zz represents the first two terms 
of the second equation in ( )68|) . On the other hand, the CP- violating contri- 
bution of the force is treated as a source S a . T^ a represents elastic scattering 
rates while Y™ ela stands for the particle number changing interactions that 
involve the chemical potentials ftp. The vectors represent which particles 
participate in a specific interaction. 

A conserved current can be represented by a vector d a . In this case all 
interactions have to preserve the current, ^ Q c a d a = 0, and the current 
should be unsourced, d a S a = 0. An example for conserved quantities are 
electric charge in the broken phase or baryon number if the weak sphaleron 
process is neglected. 

3.2 A simple diffusion network 

In order to determine the final baryon asymmetry, one has to set up a set 
of transport equations that contains all relevant degrees of freedom. The 
sphaleron rate will be one of the smallest interaction rates in the this system 
such that is suffices to neglect backreactions and determine the net baryon 
number from the left-handed particle density as described in the beginning 
of appendix |A] 

In our toy model the CP-violating source is in the top sector such that 
we first consider all fast interaction rates involving the tops. These are the 
Yukawa interactions with the Higgs, the electroweak interactions with the 
W-bosons and the strong sphaleron rate that involves all quarks. In the 
broken phase the Higgs vev induces chiral flips between left- and right-handed 
tops and also Higgs decay into W-bosons. The relevant particle changing 
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interaction rates are [3H 150] 



t L <r+t R + h 



4.2 x 1(T 3 T, 



2 

mi 



t L <-> *H 



r ss ~ 4.9 x i(r 4 T, 



t L + b L + Au L ^t R + b R + Au R 




(80) 



where and u R collectively denote the left- and right-handed light quarks. 

Next, the elastic scattering rates of the Higgs and the top have to spec- 
ified. These are usually given in terms of the diffusion constants as defined 
in (E7D and calculated in [261 E3 EEJ E3 [53] and [31] 



The Higgs and W-bosons decay quickly in the broken phase such that 
neglecting their chemical potential does not have a large impact on the final 
baryon asymmetry. A detailed analysis concerning this point can be found 
in [49] . Furthermore, the interactions with the W-bosons are rather fast such 
that left-handed up and down quarks have similar chemical potentials. The 
right-handed bottom and the light quarks are only sourced by the strong 
sphaleron rate and otherwise interact only with very small Yukawa interac- 
tions. Hence, the chemical potential of the light right-handed quarks equals 
the one of the right-handed bottom quark while the light left-handed quarks 
have the opposite chemical potential. 

Up to this point, the remaining degrees of freedom are the left-handed 
top and bottom quark with chemical potential /x 9 , and the right-handed top 
and bottom quarks denoted fi t and respectively. The light right-handed 
quarks have the same chemical potential as the right-handed bottom quark, 

and the light left-handed quarks the opposite chemical potential. 

Conservation of baryon number then relates these chemical potentials as 



The light quarks cancel in this equation since left- and right-handed particles 
have opposite chemical potentials. We also neglect the bottom masses, K h = 
kq = 1. The chemical potential of the right-handed bottom can then be 




(81) 



(fit + l)n q + fl b + KtUt = • 



(82) 
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eliminated in the remaining network. For example, the strong sphaleron 
couples to the combination 

2fi q -fi t - 9fi b = (9«t + 11K + (9«* - l)m . (83) 

In the first line, the term — 9/i& represents the nine light quark chiralities in- 
cluding the right-handed bottom. Ultimately, the left-handed baryon chem- 
ical potential entering the sphaleron process (see eq. ( \122\ of appendix [A} is 
given by 

A*l = 2/ib = (3 + 2Kt)(jL g + 2KtfH ■ (84) 

The contribution —2^ represents the left-handed quarks of the two light 
families. Notice that if the top is assumed to be light, n t = 1, the combination 
of the chemical potentials that enters the weak sphaleron process ( 183]) is 
proportional to the combination of chemical potentials that enters the strong 
sphaleron process f)84p . Hence, in this limit the final baryon asymmetry is 
suppressed by the strong sphaleron rate [53] . 

We do not quote the full set of equations here. The explicit equations for a 
network including the Higgs and W-boson fields can e.g. be found in refs [49J. 
The reduced network without Higgs field has been used in ref. [55] and [56] . 
A generalization to the two Higgs doublet model is given in ref. [57]. Also the 
generalization to supersymmetric extensions is extensively discussed in the 
literature. This includes new damping rates [58] but also much more compli- 
cated diffusion networks. In many cases it is assumed that super-gauge inter- 
actions are in equilibrium such that particle species and their superpartners 
share the same chemical potential. If this assumption is relaxed, the outcome 
of the diffusion network depends on many more parameters as e.g. the mass 
spectrum of all the superpartners. This can lead to very large correction and 
even to a change in sign in the final baryon asymmetry [59] 1601 |6T] . 

In the following, we present some results from [56]. In order to provide 
the results as model-independent as possible, the source in the top sector has 
been parametrized via the mass term as 

m t = y t (f>{z) e @t(z) , 

using 

Hz) = ^(l + tanh(z/C)) , 
S t (z) = -^(l+tanh(z/C)) • (85) 
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The final asymmetry is then proportional to the change in the top mass 
phase during the phase transition, AQ t - Otherwise, it only depends on the 
dimensionless quantities C /T C and £ W T C . 

Figure [3] shows the required change in the top mass phase during the phase 
transition AQ t in order to reproduce the observed baryon asymmetry. The 
baryon asymmetry is very sensitive to the strength of the phase transition, 
4> C /T C . Furthermore, as expected a larger wall thickness reduces the produced 
asymmetry. For phase transitions that barely fulfill the baryon washout 
criterion, C ~ T c , a change of phase of order AQ t ^ 0.3 — 0.6 is required for 
realistic wall thicknesses, £ w T c ~ 2 — 8. 
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Figure 3: The plots show the required change in the top mass phase during the 
phase transition AQt in order to reproduce the observed baryon asymmetry. In 
the upper plot the wall thickness in terms of the temperature is kept constant, 
while in the bottom plot the wall thickness in terms of the critical vev is kept 
constant. The plots are adapted from ref. [56] . 
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4 Models 



The crucial ingredients of electroweak baryogenesis are a strongly first-order 
phase transition and an appropriate source of CP violation. 

A strong electroweak phase transition is needed for several reasons. First, 
the nucleated bubbles during the first-order phase transition are the source 
that drives the plasma locally out of equilibrium and facilitates the establish- 
ing of sizable CP-violating currents. Second, the baryon number violating 
sphaleron processes have to be sufficiently suppressed after the phase transi- 
tion in order to avoid the washout of the just produced baryon asymmetry. 
This leads to a constraint on the Higgs vev C and the phase transition tem- 
perature T c of the form (see sec. |A] for a short derivation of this bound) 



In the Standard Model, it is well known that a first-order phase transition 
is excluded for Higgs masses beyond roughly the W-boson mass [62]. To 
fulfill the washout criterion ( )86|) . a Higgs mass below ~ 40 GeV would even 
be necessary. This is in contrast to the Higgs mass bound from LEP of 
mil > 114 GeV. Generally, a strong phase transition fulfilling ( 186|) requires 
either an extended scalar sector or at least new degrees of freedom that are 
strongly coupled to the Higgs. 

In electroweak baryogenesis, an appropriate source of CP violation has 
to be in the form of a complex mass matrix that changes during the phase 
transition such that a chiral flux is generated close to the bubble wall. This is 
achieved by coupling the corresponding particles to a vev that constitutes the 
nucleating bubbles of the phase transition. In many cases this vev arises from 
the physical Higgs field but more complicated scalar sectors tend to increase 
the prospects of electroweak baryogenesis. This is because the masses of the 
Standard Model fermions are proportional to the Higgs vev. Comparison 
with the sources in (148]) then shows that CP violation is absent. Hence 
either the masses of the Standard model fermions need to be modified or a 
new fermionic particle is responsible for the CP-violating flux. In the latter 
case, the CP-violating flux has to ultimately bias the sphaleron rate. Hence 
this new degree of freedom is in many models charged under SU(2)l- 

At the same time, these new features can leave traces in collider and low 
energy probes. One major constraint comes hereby from electric dipole mo- 
ments that constrain new sources of CP violation. Often the induced electric 




(86) 
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dipole moments arise only at two loop. Still, current bounds on the electron 
EDM, d e < 1-05 x l(r 27 ecm [63], and neutron EDM, d n < 2.9 x 10" 26 ecm 
[64] , heavily constrain realistic models of electroweak baryogenesis. Also the 
new degrees of freedom responsible for a strong phase transition can have 
measurable implications. The prime example for this is the Minimal Su- 
persymmetric Standard Model (MSSM) where only light right-handed stops 
can yield a sufficiently strong phase transition. Such light stops would be 
copiously produced at the LHC what leads to additional constraints. 

In the following we discuss several models in which electroweak baryogene- 
sis is feasible. We start with relatively simple models with higher dimensional 
operators and the two Higgs doublet model in which the semi-classical force 
is operative. Then we discuss the MSSM and its extensions that requires a 
more sophisticated treatment of CP violation from flavor mixing. 

4.1 The Standard Model with a low cutoff 

From a bottom-up perspective, the minimal approach to extensions of the 
Standard Model is to insist on the particle content of the Standard Model 
and only extend the Lagrangian by higher dimensional operators. Since elec- 
troweak baryogenesis requires sizable deviations from the Standard Model at 
around the weak scale, the suppression of the higher-dimensional operators 
and the physical cutoff of the theory cannot be much larger in this frame 
work. Still, the higher dimensional operators can have an important im- 
pact on the phase transition, provide new sources of CP violation and make 
electroweak baryogenesis a viable option. 

4.1.1 Phase transition 

The leading operator that modifies the Higgs potential is of the form (<3>t$) 3 , 
such that the scalar potential of the Higgs vev reads 

V{<j>) = /iV + A0 4 + . (87) 

The new scale A is the cutoff of the theory where new degrees of freedom 
become relevant or at least strong coupling phenomena occur. This form 
of potential can lead to a strong phase transition already in the mean-field 
approximation where temperature effects only contribute to the quadratic 
Higgs term, AVr — cT 2 cf) 2 . The barrier is then produced by balancing a 
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negative quartic term, A < 0, with the positive 6 operator [651 EE]. The 
critical temperature is then 

_ A*m% + 2A 2 m^ - 30* 
Tc " 16^ ' (88) 

where the parameters \i and A have been expressed in terms of the physical 
Higgs mass m# and the observed Higgs vev 0o — 246 GeV. The critical vev 
is given by 

* = < 89) 

There is also an upper limit on A where the phase transition becomes second 
order and a lower bound from the fact that the broken phase is the global 
minimum at T = 0. As usual an increase in the Higgs mass makes the phase 
transition weaker. The washout criterion flBBl) translates into an upper bound 
on A. In the full one-loop analysis, the values are [66] 

A < 800 GeV, m H = 125 GeV , 

A < 900 GeV, m H = 115 GeV . (90) 

A peculiar feature of the model seems to be that the coefficient of the quar- 
tic A is negative. However, a negative quartic can arise quite naturally in 
effective actions, for example when a heavy scalar is integrated out [65J. 

4.1.2 Electroweak baryogenesis 

Electroweak baryogenesis was considered for this model in [55]. As an addi- 
tional efficient source of CP violation served a dimension-six coupling between 
the Higgs $ and the up-quarks 

C3 -^ L (&®)q i ®u j + h.c, (91) 

^CP 

in combination with the usual Yukawa coupling 

£ 3 yijqi $ uj + h.c. , (92) 
The resulting fermion masses during the phase transition read 
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Figure 4: The wall thickness £ w as a function of the Higgs mass. The plot shows 
also the corresponding values of the scale of new physics A and the ratio £ = 4> C /T C . 
Plot adapted from [55] . 

what leads to a CP- violating semi-classical force if there are relative complex 
phases between and x^-. The most important effect is in the top sector, 
since the other quarks are too light to yield a sizable CP-violating flux along 
the bubble wall. The change of the phase is hence of order 



where x t denotes the 33 element of the Xij coupling in the mass eigenbasis of 
the quarks. 

The system of transport equations is the one discussed in section I3.2I 
The only degrees of freedom are the ones from the Standard Model and the 
dominant source of CP violation is the semi-classical force in the top sector. 
The only missing ingredient in the present context is the wall thickness. A 
numerical analysis of all the characteristics of the phase transition (see Fig. H]) 
and the analysis of the produced baryon asymmetry is given in ref. [55]. The 
final baryon asymmetry is very sensitive to the scale A. The main influence 




CP 



(94) 
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comes hereby from the relation between A and the critical Higgs vev <p c . The 
semi-classical force (IBTl) is proportional to 4>l via the top mass and another 
factor (fi 2 c stems from the change in the phase (1941) . Besides, the wall thickness 
£ w tends to be smaller for stronger phase transitions and hence lower values 
of A. For ^s(xt) 1, and imposing the observed baryon asymmetry implies 
the bound A CP ~ A < 650 GeV. 



4.1.3 Collider and low energy probes of the model 

Since the model does not contain any new degrees of freedom, no spectacular 
signatures are expected at colliders. Still, the higher dimensional operators 
can lead to measurable deviations from the Standard Model. 

In connection to the phase transition, the new operator 6 is the essential 
ingredient. The main collider trace of this new operator is a deviation of the 
self-couplings of the Higgs in terms of the Higgs mass [55]. The deviations 
from the Standard model couplings read 

Ax = 3— + 6-, ^ = 3^ + 36-, (95) 

where \i (rj) denote the cubic (quartic) self-coupling of the Higgs field. The 
deviations are pronounced for small Higgs mass, e.g. \i ~ 2/i 5M for m# = 125 
GeV and A = 650 GeV. Still, the discovery of a deviation of this size requires 
a linear collider [65]. However, in combination with EDM bounds, viable 
baryogenesis requires an even stronger phase transition what makes even 
larger deviations in the Higgs sector necessary. This is discussed next. 

The new source of CP violation gives potentially much stronger bounds 
in light of observed limits on flavor changing neutral currents. However, 
these bounds are more model dependent and in particular hinge on the fla- 
vor structure X{j of the new operator (l9Tj) . Flavor changing neutral currents 
potentially arise, because the mass term f )93|) is not proportional to the cou- 
pling between the Higgs and the fermions 



1 3v 2 
Yij = Vij V2 +Xij VSA^' (96) 



If the couplings were random numbers of order unity, large deviations in 
the first two quark families could be observed. For example, the operator 
f )9TT) would affect K — K mixing [55] what implies a bound Aqp > 10 7 GeV. 
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Figure 5: Two loop contribution to the electron EDM of Barr-Zee type. 

If on the other hand, Xij has a similar flavor structure as jjij, the model is 
consistent with these constraints as long as Aqp > 500 GeV. Such a setting 
is well motivated from the hypothesis of minimal flavor violation [67J and 
can be achieved in Froggatt-Nielsen type models. 

Constraints from the electric dipole moments are more robust, since they 
occur even in a one flavor scheme with a relative phase between the couplings 
yt and Xt- Since this phase is also essential for the production of the baryon 
asymmetry, this provides a direct link between low energy observables and 
baryogenesis. The dominant constraints [68J come from the Barr-Zee type 
contributions to the neutron and electron EDMs (see Fig. [5]). For a Higgs 
mass of m H = 125 GeV this constraint reads A C p > ^/Q(x t ) X 750 GeV and 
becomes slightly weaker for larger Higgs masses. 

4.1.4 Summary 

Electroweak baryogenesis is a viable possibility in the Standard Model with 
a low cutoff. The strongest constraints on the model come from the cubic 
Higgs self-coupling and the upper bounds on the neutron EDM. On general 
grounds one expects that the new operator in fl57|) that makes the phase 
transition strong and the operators that provide the CP violation (|9~T|) are of 
similar size, A ~ Aqp- This is indeed possible for A somewhat smaller than 
Ac p. A possible set of parameters is for example 

A ~ 500 GeV, A CP ~ 1000 GeV, $(x t )~l. (97) 
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If the Higgs is rather light (m# ~ 125 GeV), this particular set of parameters 
will be tested in the near future. The cubic coupling is enhanced by a factor 
~ 3 compared to the Standard Model what could even be in reach for the 
high luminosity LHC Furthermore, the next generation of neutron EDM 
measurements (assuming an improvement of factor 10 in sensitivity) can 
rule out this model of electroweak baryogenesis. 

4.2 Low cutoff: Singlet extension 

The best motivation for extensions of the Standard Model with a low cut- 
off comes from composite Higgs models. In composite Higgs models, the 
light spectrum of the scalar sector depends on the co-set structure of the 
strongly coupled sector. The degrees of freedom below the scale of strong 
coupling arise as bound states with pseudo-Goldstone nature. In the mini- 
mal model [69J [701 EE] the Higgs is the pseudo-Goldstone boson of the the 
breaking pattern SO (5) — >■ 50(4) where 5*0(4) represents the custodial sym- 
metry of the Higgs sector. In non-minimal models an extended scalar sector 
appears at low temperatures. In the following we discuss the model with 
the breaking pattern 50(6) — >■ 50(5) that serves as a UV completion of the 
singlet extension of the Standard Model with a low cutoff [72] . 

From a phenomenological point of view, electroweak baryogenesis can be 
more easily realized in this model than in the Standard Model with a low 
cutoff. First, the phase transition can be strong already with a renormaliz- 
able scalar potential and in the mean-field approximation and does not rely 
on higher dimensional operators at all. Second, the leading source of CP 
violation arises already at dimension five. This allows to push the cutoff to a 
few TeV what is advantageous in view of flavor physics. Last, the dominant 
contribution to EDM constraints stem from a mixing between the Higgs and 
the additional singlet degree of freedom. As long as this mixing is small, 
current constraints from low energy probes are easily fulfilled. 

4.2.1 Phase transition 

As mentioned before, the phase transition can be already strong in mean-field 
approximation with only renormalizable operators in the scalar potential. 
Interestingly, this is even true if a Z 2 -symmetry is imposed on the singlet, 
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s — s. Consider the following potential at the critical temperature: 



V 



T=T r . 



.V 



(98) 



The variables <fi and s denote the Higgs and singlet vev and <ft c and s c the 
corresponding values of the vevs in the SU{2)l and Z 2 -breaking phases at 
the critical temperature. The first term constitutes a Mexican hat potential 
with a flat direction that connects the S , f/(2)£,-breaking with the Z 2 -breaking 
phase. The second term lifts this flat direction and creates a barrier between 
the two degenerate minima of the potential. 

Thermal corrections in the mean-field approximation can be added to this 
potential via 



AV T 

where the two coefficients 



1 



{c^ 2 + c s s 2 ){T 2 -T 2 c ) 



(99) 



and c, read 



48 
1 

12 



9g 2 + 3g' 2 + 12y 2 + 24A + 4^/\\~ s + 2k 
3A s + 4 v / AA^ + 2/t 



(100) 



and we defined A s = ' s% In total the model has four free parameters that 
can be fixed using the observed Higgs vev <fi = 246 GeV, the Higgs mass, the 
singlet mass and the critical temperature. A lower bound on the singlet mass 
results from the requirement of a first-order phase transition (k > 0) while 
an upper bound on the singlet mass arises from the requirement that the 
SU (2) L-broken phase is the global minimum at T = 0. Detailed plots are 
given in [73] and also in [74]. For fixed Higgs and singlet masses, the critical 
temperature can always be reduced down to the point where the system 
becomes very strong, <p c /T c ~ a few. 

In fact, the phase transition proceeds in two stages in this model: At 
very high temperatures, the singlet vev as well as the Higgs vev vanish and 
neither the electroweak SU(2)l nor the Z 2 symmetry are broken. At lower 
temperatures the singlet develops a vev that breaks the Z 2 symmetry. De- 
pending on the parameters, this process can happen at several hundred GeV 
and is probably rather a cross over than a phase transition. At this stage, 
domain walls are generated. However, the domain walls are harmless to big 
bang nucleosynthesis since they disappear in the next stage when the system 
transits from the Z 2 -breaking phase to the electroweak breaking one. 
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4.2.2 Electroweak baryogenesis 

In contrast to the Standard Model with a low cutoff, its singlet extension 
already has a powerful source of CP violation at dimension five 

C3 _pL S q i $ Uj+ h.c, (101) 
The resulting fermion masses during the phase transition read 

(f> s 4> 

m *3 = Vij ^= + By ^2 kcp ' ( 102 ) 

what again leads to a CP-violating semi-classical force if there are relative 
complex phases between and x^. Following the rationale of the Standard 
Model with low cutoff, we focus on the top sector. The change of the phase 
of the top mass is of order 

AQ t ~ 9f(z t ) , (103) 

where x t denotes again the coupling in the mass eigenbasis of the quarks. 
Compared to the minimal model with cutoff, the singlet extension has several 
nice features in view of baryogenesis. First, the phase transition can be rather 
strong without coming into conflict with a low cutoff. Next, the change of 
phase (11031) is only suppressed by one power of A what makes baryogenesis 
in this model easier compatible with a cutoff A ~ 2 — 3 TeV. With such a 
high cutoff, it is e.g. possible to solve the flavor problem using the 5D GIM 
mechanism in specific realizations of the composite Higgs mechanism [75j. 
Furthermore, the singlet vev is in principle expected^ to be larger than the 
Higgs what further increases the source (11031) . Some numerical results are 
shown in Fig. |6j Electroweak baryogenesis can be viable for AQ t > 1 what 
translates into the bound A^p < a few TeV. 

Notice that if the scalar potential is completely Z 2 symmetric the baryon 
asymmetry is suppressed. As mentioned above, domain walls are generated 
at intermediate scales where the singlet vev breaks the Z 2 spontaneously. At 
this stage the Universe is divided into regions with positive/negative singlet 
vev. These regions produce opposite baryon numbers during the electroweak 
phase transition. In order to avoid this problem, the Z 2 has to be slightly 
broken. Already a very small breaking leads to a disappearance of the domain 
walls and preserves the baryon asymmetry |56| . 



7 5D realizations of the composite Higgs require a slight tuning to make the electroweak 
scale and hence the Higgs vev small |70j . 
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Figure 6: The shaded region shows possible models in the plane £ = (j) c /T c 
versus the change in top phase AQf The Higgs and singlet masses are ran = 120 
GeV and m s = 130 GeV. The scale of new physics is Aqp = 500 GeV. The red 
lines denote the parameters that reproduce the observed baryon asymmetry. Plot 
adapted from |56j . 

4.2.3 Collider and low energy probes of the model 

Unlike the Standard Model case, the additional CP-violating operator does 
not give rise to dangerous flavor observables. First, if the model is approx- 
imately Z 2 -symmetric, the operator (11 01 1) is absent after the electroweak 
phase transition. Even if the scalar field has a (small) vev after the elec- 
troweak phase transition, the Yukawa interactions with the fermions can be 
diagonalized simultaneously with the fermionic mass terms fll02[) what sup- 
presses flavor changing neutral currents to higher loop order. 

In terms of collider traces and electric dipole moments, deviations from 
the Standard Model arise mostly from a singlet-Higgs mixing. As mentioned 
above, a very small Z 2 breaking is required for viable baryogenesis but it is 
easily compatible with bounds from electroweak precision tests or EDMs as 
seen in Fig. [7J 

Another characteristic signal of the model would be a Higgs decay into 
four fermions via two singlets. If this process can be tested depends however 
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m / GeV 



Figure 7: Bounds on the mass rri2 of the mostly-singlet mass eigenstate, and the 
mixing angle 9. The mostly-Higgs state has a mass mi = 120 GeV. The scale of 
new physics is Acp = 500 GeV. Plot adapted from [56] . 

on the coupling of the singlet to the (non-top) fermions and on the mass of 
the singlet. 

4.2.4 Summary 

Electroweak baryogenesis can be very easily realized in singlet extensions of 
the Standard Model with a low cutoff. Since efficient sources of CP violation 
are present with dimension five operators, the cutoff scale can be slightly 
larger than in the case of its minimal cousin, A < a few TeV. Also collider 
bounds and low energy probes can be easily avoided if the Z 2 symmetry of 
the singlet sector is only weakly broken. This makes the model compati- 
ble with phenomenology and insensitive to EDM constraints. However, one 
has to notice that this is somewhat against the philosophy of electroweak 
baryogenesis that has falsifiability at its core. 
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4.3 Two-Higgs-doublet model 

In the two Higgs doublet (THD) model, all necessary ingredients of elec- 
troweak baryogenesis are present, even if only renormalizable operators are 
considered. The most general scalar potential reads 



The potential contains two complex (potentially CP-violating) couplings 
/i 3 e* Q and A5. Following the conventions of [57J, we choose A5 to be real 
such that a parametrizes CP violation in the scalar sector. As we will see in 
the next section, the complexity of the scalar potential is also high enough 
to provide a strong first-order phase transition. 

4.3.1 Phase transition 

In principle there are two regimes in parameter space with a strong first- 
order phase transition. The first one is similar to the case discussed in the 
singlet extension of section 14.21 The phase transition again proceeds in two 
steps, but unlike in the singlet extension, already this first phase transition 
breaks the electroweak symmetry in the THD model. This implies that for 
viable electroweak baryogenesis, this first phase transition has to be strongly 
first-order, which is not so easily achieved. We hence dismiss this possibility 
of a two-stage phase transition in the following. 

The reason that the phase transition can be much stronger than in the 
Standard Model is two-fold. The first is that both Higgs doublets acquire a 
vev after the phase transition and the form of the potential implies that the 
ratio tan /3 of these two vevs 



<$i) = (0, he 101 ), ($ 2 > = (0, h 2 e i02 ), tan/3 = V^2, (105) 



is not constant during the phase transition. The potential in terms of the 
vev 2 = h\ + h\ is hence not necessarily polynomial and eventually develops 
a barrier between the two minima at the critical temperature. The second 
reason is that the scalar potential has enough free parameters to decouple 




(104) 
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Figure 8: Lines of constant £ = (f> c /T c and t yj CIS 3j function of the two scalar 
masses and ran and for ^ = 100 GeV and a = 0. Plot adapted from [57] . 

the Higgs mass from the quartic coupling, that in the Standard Model are 
related by m 2 h = 2A0q. It is hence possible to obtain a strong phase transition 
from the thermal cubic contributions to the effective potential and to satisfy 
at the same time the LEP bound of rrih > 114 GeV. 

Overall, relatively strong phase transitions, £ = 4> C /T C > 1.5, i w T c < 10, 
are possible for a Higgs mass above the LEP bound [761 EH EH EH]- Some 
examples are shown in Fig. El 

4.3.2 Electroweak baryogenesis 

The most general THD model with Yukawa couplings of fermions to all two 
Higgs fields suffers from flavor changing neutral currents already on tree level. 
To avoid this problem, usually an additional Z 2 symmetry is invoked that 
allows to couple the fermions only to one of the two doublets 



$i -> -$i , d -> ±d, 



(106) 
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where depending on the sign in the down sector the THD models type I and 
type II result, respectively. Notice that the complex phase a in the potential 
(11041) breaks this symmetry explicitly such that electroweak baryogenesis is 
not possible if this symmetry is also imposed on the scalar sector. 

As before, the main source of baryogenesis comes from the top sector and 
the corresponding Yukawa coupling is of the form 

C3y t Q^ 2 t. (107) 

The CP- violating source comes in this model not from the interplay between 
two operators that both contribute to the top mass, but from the change of 
the complex phase 2 in the Higgs field that couples to the top 

m t = ^Lh 2 e ie2 . (108) 
v2 

The change of 02 during the phase transition is hereby induced by the de- 
pendence of the scalar potential on the relative phase AO = (©i — 02)/2 
that arises in the contributions involving a. 

In [57] a part of the parameter space of the THD model is analyzed 
under the assumption that tan f3 does not change during the phase transition. 
However, using this assumption can lead to over-estimating the present CP 
violation as detailed in [79]. The reason is the following: The kinetic terms 
of the Higgs fields coming from the phases yields in the effective action for 
the vevs the contributions 



s 3 kQ[) 2 hi + ke' 2 ) 2 hl 



= I(e') 2 (^ + ^) + i(A6') 2 (^ + ^) + Ae'e'(^-^), (109) 

where we defined the average phase G = (0i + 02)/2 and the relative phase 
A6 = (6i — 6 2 )/2. Since the effective potential does not depend on the 
average phase 9, one finds (using the equations of motion) 

e ' - "frf Ae ' • < 110 » 

Reinserting this into the kinetic term gives 
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and for the individual phases 

2h 2 -2h 2 



ei = if^Ae', e' 2 = T ^Ae'. (112) 



From this it follows that CP violation in the top sector vanishes if one of 
the vevs vanishes in the trajectory during the phase transition (for hi = 
the top mass vanishes while for h\ = the phase B2 is constant and no 
semi-classical force is present). On the other hand, one can always make a 
basis choice where only one of the Higgs fields has a vev in the broken phase. 
If a constant tan/3 was imposed in this basis, CP violation would be absent. 
So the assumption of constant tan (3 is not only a basis dependent statement 
but also of major importance for CP violation. Furthermore, this argument 
shows that the baryon asymmetry should be suppressed in the limit of very 
large or very small tan f3. 

Numerically, the study [57] found that a baryon asymmetry a few times 
larger than the observed one is possible in this setup. In contrast, the analysis 
[79] additionally implemented (very strict) constraints on Z — > bb and found 
generically a smaller baryon asymmetry. 



4.3.3 Collider and low energy probes of the model 

The THD model and its collider phenomenology is widely studied in the 
literature (for a recent review see [80J ) . In the context of electroweak baryo- 
genesis the main signatures are again the electron and neutron EDMs but 
also the masses of the additional Higgses that have a large impact on the 
strength of the electroweak phase transition. 

The study [75] found in agreement with [57] that for fixed Higgs mass 
rrih stronger phase transitions can be obtained especially if the additional 
Higgses are rather heavy. As explained in [57] , this arises from the fact that 
the larger masses stem from larger quartic couplings and hence corresponds 
not to a decoupling of the additional Higgses. On the other hand, the quartic 
couplings are not so essential for collider searches and EDM constraints such 
that in this limit electroweak baryogenesis is rather unconstrained in the 
THD model. The limiting factor in this regime is that one wants to preserve 
perturbativity of the quartic couplings. 

As mentioned before, additional constraints come from Z — > bb. The 
main deviation from the Standard Model stems from the loop contributions 
of the charged Higgses to this process. In general, this drives the model to 
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larger masses of the charged Higgses and to larger tan /3. This is problematic 
for electroweak baryogenesis, since large values of tan/3 suppress the CP- 
violating semi-classical force. In [79] very strict bounds on this process (i.e. 
66% C.L.) have been implemented what has a large impact on the final 
baryon asymmetry. If this constraint is treated more permissively (e.g. with 
95% C.L.) the corresponding bound is not so relevant and tan (3 is relatively 
unconstrained. 

4.3.4 Summary 

Electroweak baryogenesis is a viable option in the THD model. Without tun- 
ing the model allows for a strong first-order phase transition and sufficient 
CP violation in the scalar sector consistent with EDMs and collider probes. 
The main disadvantage of the model is that it does not have many benefits 
beyond electroweak baryogenesis. In particular, the hierarchy problem re- 
mains unsolved and flavor issues cannot be solved by a discrete symmetry in 
the cases where electroweak baryogenesis is possible. 

Over all, an improvement of the measurement of the neutron EDM by a 
factor around ten can exclude electroweak baryogenesis in the THD model. 

4.4 MSSM 

The minimal supersymmetric standard model (MSSM) is one of the most 
widely studied models today and one of the biggest contenders for the ques- 
tion how the large hierarchy between the electroweak and the Planck scale 
can be explained. 

The analysis of electroweak baryogenesis is in the MSSM very different 
compared to other models. First of all, there is no CP violation in the 
scalar potential and the top sector (beyond the CKM CP violation of the 
SM), such the CP violation has to arise from a different source than in the 
cases discussed so far. In addition, it is not easy to obtain a strongly first- 
order phase transition in this setup. In particular, the ratio C /T C even in 
most optimistic scenarios barely fulfills the washout bound (1861) and the wall 
thickness is rather large, £ W T C ~ 20. This leads to a situation where the 
semi-classical force falls short to explain the observed baryon asymmetry. 
Hence electroweak baryogenesis in the MSSM has to be based on a different 
source of CP violation as e.g. the mixing between different charginos (and 
eventually neutralinos) that can be resonantly enhanced. 
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A more extensive recent review of electroweak baryogenesis in the MSSM 
is given in Ref. [81J and we just present a short overview of the main points 
here. 



4.4.1 Phase transition 

The scalar potential in the MSSM is much more constrained than the one of 
the general THD model. On tree level it reads 

2 I 12 

V = m\h\ + m\h\ + 2m\h x h 2 + - y {h\ - h 2 2 ) . (113) 

With this potential the mass of the lightest Higgs bosons is constrained to 
be below the Z-boson mass. This is not compatible with the bounds from 
LEP and calls for large one-loop contributions to the Higgs mass 
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The dominant contributions to the Higgs mass come hereby from the tops and 
stops that have Yukawa couplings of order one and the masses, m t = yth 2 , 

* \y t (A t h 2 - fih) ml + yfhl J v ' 

where mu, rriQ and A t are soft supersymmetry breaking terms and \x stems 
from a term in the superpotential of form W 3 fi H\ ■ H 2 . In order to obtain 
a Higgs mass m h ~ 125 GeV, at least one of the stops has to be rather heavy, 
mi > 30 TeV. This can be achieved by either a large soft mass m,Q or by a 
large off-diagonal contribution from the A t term. 

The second option is not compatible with a strongly first-order phase 
transition as we will see in the following. As in the standard model, the po- 
tential barrier that is responsible for the first-order phase transition can only 
arise from thermal cubic terms in the effective potential (see appendix ID]) . 
Besides the degrees of the freedom of the Standard Model, only the stops can 
give such a sizable cubic term [821 [83j EU |85l |86] . This means in turn that 
the right-handed stop (that is less constrained by electroweak precision tests 
than its left-handed partner) has to be very light. In particular, a cubic term 
is only delivered if the mixing between the stops is small and the thermal 
mass of the right-handed stop is countered by a negative soft mass, i.e. 



IJT) = mfj + yX + n(T) iR ~ yX . (116) 
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Figure 9: The window of strong enough phase transition, 4> C /T C > 1.0, in the 
Higgs mass versus light stop mass plane for the MSSM. A strong phase transition 
and a Higgs mass ~ 125 GeV can only be achieved at the cost of a very heavy 
left-handed stop, tuq ~ 10 6 TeV. Plot adapted from [87] , 

Additional constraints arise from the requirement that tan /3 is not too large 
and that the stop do not develop a vev at low temperature what would lead 
to a spontaneous breaking of color. The results of this analysis from [87] 
is shown in Fig. HI These results also have been qualitatively confirmed in 
lattice calculations [88J. 

4.4.2 Electroweak baryogenesis 

As alluded in section 12.51 the determination of the baryon asymmetry in the 
MSSM is a controversial topic. One difference to the other models discussed 
so far is that CP violation does not arise in the top sector. The dominant 
source of CP violation turns out to be the charginos and neutralinos. For 
example the chargino mass can be written 



(117) 
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where M2 and /x can contain a complex phase. 

This mass matrix will lead to a source of the semi-classical force type 
according to eq. fl50|) . However, the phase transition in the MSSM is rel- 
atively weak [891 [90], C /T C ~ 1, £ W T C ~ 20, such that this source of CP 
violation is not sufficient to explain the observed baryon asymmetry once 
EDM constraints are imposed. 

Hence baryogenesis has to be driven by mixing effects in the MSSM. Para- 
metrically, mixing effects are less suppressed because they appear already at 
first order in gradients as seen in eq. (150]) . The determination of the baryon 
asymmetry based on these mixing effects is to certain extent still an open 
issue. The mass insertion formalism yields very large baryon asymmetry [UJ 
but suffers from conceptual problems (see sec. 12. 6p . Part of these problems 
can be overcome by resumming Higgs insertions [42] but also in this frame- 
work some issues concerning finiteness of the results and how transport is 
established remains. Conceptually the cleanest way to tackle this problem is 
to use the first principle approach in the Kadanoff-Baym framework. This 
was done in the analysis [21] that particularly highlighted the importance of 
flavor oscillations. But also in this study many simplifying assumptions have 
been used. Namely, the coherent off-diagonal densities have been assumed 
to be small. In particular, all contributions that are nominally second order 
in gradients have been neglected. If these contributions are really small is 
not so clear since resonant effects can become important when the oscillation 
length is close to the wall thickness [20]. Naively, this resonance condition is 
for the MSSM charginos only fulfilled for rather hard modes (that are sparse 
in the plasma) but this does not guarantee that the resonance can give a 
large enhancement of the baryon asymmetry. 

But there are also some features that are shared by all approaches. For 
example the baryon asymmetry is suppressed when the charginos are not 
almost mass degenerate or have a mass much larger than the temperature. 
This is seen in Fig. [TO] that shows the regions of viable baryogenesis as a 
function of the two chargino mass parameters. A selection of quantitative 
results of chargino driven baryogenesis in the MSSM is collected in Table [TJ 

Beyond theses studies, neutralino [HI] or stop driven [92] baryogenesis was 
considered for the MSSM in the literature. Neutralinos have the advantage 
that they do not suffer from as large EDM constraints as charginos but also 
are somewhat less efficient in producing the baryon asymmetry [91J. 



4 MODELS 



54 



m A = 2 00 GeV m A = 400 GeV 




100 200 300 400 100 200 300 400 

M2 / GeV M2 / GeV 



Figure 10: Contours of the regions with viable baryogenesis as a function of 
the two chargino mass parameters M2 and \x. In the black region the baryon 
asymmetry is larger than observed. Plot adapted from |24j . 



paper 


method 


r)/ Vobs 


gU (2000) 


mass insertion formalism; no Higgs re- 


~ 35 




summation 




pE2] (2002) 


mass insertion formalism; including 


~ 10 




Higgs resummation 




H3] (2004) 


mass insertion formalism; no Higgs 


~ 140 




resummation; more realistic diffusion 






network 




[21] (2005) 


Kadanoff-Baym formalism; flavor oscil- 


~ 3.5 




lations; assumes the adiabatic regime 





Table 1: The largest possible baryon asymmetry for almost mass degenerate 
charginos and a maximal CP-violating phase. 

4.4.3 Collider and low energy probes of the model 

In the context of electroweak baryogenesis, the MSSM provides some special 
signatures. The first class of signals comes from the new source of CP viola- 
tion in the chargino sector. Since the charginos cannot be much heavier than 
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Figure 11: Contours of the electron and neutron EDMs as functions of the two 
chargino mass parameters and for a maximal CP phase. Plot adapted from |93| . 

the electroweak scale in electroweak baryogenesis, this leads to Barr-Zee type 
contributions to the neutron and electron EDMs that are sizable and can be 
already in conflict with experimental bounds. Furthermore, the dependence 
of the electron EDM on tan/3 and the chargino masses is quite similar to 
the dependence of the baryon asymmetry j93j Ell [95] (see Fig. [TT]) . This 
implies that the complex phase in the chargino sector cannot be larger than 
arg(yU*M 2 ) < 0.05. This excludes chargino driven electroweak baryogenesis 
in the MSSM in the most conservative approaches (see Tabled]). 

The second class of constraints is connected to the requirement of a strong 
first-order phase transition. The most severe is hereby the occurrence of stops 
close to the LEP bound [96] |97] [98] - More recently, direct searches at LHC 
are sensitive to light stops such that this is only viable if stop decays are 
concealed through neutralino states with similar mass or some alternative 
mechanism [87J. Still, the light stops would have a large impact on the 
Higgs search. In particular, they increase the loop-induced Higgs production 
rate by gluon fusion by a factor 2 to 3. Besides, light stops lead to a reduced 
branching ratio for Higgs to di-photons due to a destructive interference with 
the dominant W-boson loop. Overall, light stops lead to an enhancement of 
the rate gg —> H — >■ VV and a slight reduction of the rate gg — > H — > 77 
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compared to the Standard Model. This produces a tension with the current 
data from Higgs searches (99J 11001 HOlj that can be partially relaxed by 
further assumptions about the particle spectrum [87J. 

4.4.4 Summary 

Electroweak baryogenesis in the MSSM is an appealing scenario because the 
MSSM is the minimal setup that solves the hierarchy problem in a perturba- 
tive framework. A Higgs mass of rrih ~ 125 GeV produces a tension with min- 
imal supersymmetric models, particularly when a strongly first-order phase 
transition is demanded. This requires in addition very light stops right above 
the LEP bound. That these states have been missed at LHC so far is possible 
but only in case of peculiar masses for the particles that appear in the decay 
chain of the stops [ST]. Also the EDM constraints are generically in conflict 
with chargino driven electroweak baryogenesis. So either other sources of 
CP violation (e.g. neutralinos) have to be utilized or the EDMs are small 
because of a cancellation of different contributions. 

In summary, there remains a region of parameter space in the MSSM in 
which electroweak baryogenesis is still viable. Still, this possibility appears 
rather contrived with several requirements arising in different sectors. How- 
ever, the most constraining requirements can be traced back to the fact that 
Higgs masses of ~ 125 GeV are not easily realized in the MSSM. In ex- 
tensions of the MSSM where the Higgs mass is achieved more naturally, also 
the prospects of electroweak baryogenesis are much better. This is explicitly 
seen in the next model. 

4.5 Next-to-MSSM 

The main aim of singlet extensions of the MSSM is two-fold. First, the /i- 
problem of the MSSM is solved. This is accomplished by adding a term 
of form XSHi ■ H 2 to the superpotential. When the singlet acquires a vev 
by spontaneous symmetry breaking, this operator produces an effective /i 
term. Second, additional contributions to the lightest Higgs mass improve 
the consistency with current collider constraints. In the following we discuss 
a variant with only trilinear coupling to the Higgses and a linear term for 
the singlet in the superpotential as done in |102[ I103[ 1104] . More general 
models can lead to new phenomena as e.g. transitional CP violation [105] or 
interesting dark matter phenomenology [106] . 
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4.5.1 Phase transition 

In this model the scalar potential reads 

V = mih* + mlhl + 2mlh 1 h 2 + ^^(hl-h 2 2 ) 2 

8 

+ m 2 s s 2 + ^hjhj + a x s h\h 2 + t 8 s e i6s + h.c. . (118) 

where we defined the vev of the scalar field as (S) = s e l&a / \/2- Here, the 
parameter A results from the term \SH\ ■ H 2 in the superpotential and t s and 
a a are soft SUSY-breaking terms. Of special importance is the contribution 
Xhjhl/A which lifts the D-flat direction of the MSSM and can give a sizable 
contribution to the lightest Higgs mass. 

The phase transition can become strong due to the interplay of the singlet 
and the Higgs vevs and does not rely on thermal loop corrections. Already 
on tree level the model develops a first-order phase transition when [103J 
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(119) 



where we defined 



A 2 = ^ sin 2 (2/3) + g2+ o 9 ' 2 cos 2 (2/3) . (120) 
4 8 

For moderate values of A, Higgs masses of order nth ~ 125 GeV are possible 
and consistent with a strong phase transition. However, the parameter A 
eventually develops a Landau pole at not too high scales what implies the 
rough bound A < 0.7. 



4.5.2 Electroweak baryogenesis 

Electroweak baryogenesis is easier to realize in the nMSSM than in the MSSM 
for several reasons. First of all, the phase transition can be much stronger. 
This gives a considerable enhancement in the CP- violating source, that is 
very sensitive to <f> c /T c , but also due to a reduced wall thickness. Further- 
more, additional complex phases in the parameters t s and a\ lead to new 
sources of CP violation. In particular, the phases of the singlet and the 
Higgs fields changes during the phase transition |104[ 1107] . The former leads 
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to an additional semi-classical source in the chargino sector via the modified 
chargino mass matrix 



but also to a source in the top sector due to a change in 62 during the phase 
transition. These contributions arise in the semi-classical force approach and 
do not rely on mixing. Additional sources by mixing can be as large as in 
the MSSM but since the semi-classical forces do not require almost mass 
degenerate charginos these contributions are typically much smaller. This 
allows for a rather reliable determination of the baryon asymmetry compared 
to the MSSM. 

4.5.3 Collider and low energy probes of the model 

Compared to the MSSM, collider and EDM constraints are easier to fulfill 
in the nMSSM. As mentioned before, the lightest Higgs can obtain sizable 
mass contributions from the coupling to the singlet. However, Higgs masses 
of rrih ~ 125 GeV that rely solely on this coupling lead to a Landau pole 
in the coupling A below the GUT scale. Hence loop corrections from the 
stops and tops still have to be sizable and stops heavier than a TeV are 
required. Notice that light right-handed stops are not essential for a first- 
order phase transition, such that they can have masses similar to their left- 
handed counterparts. 

Constraints from EDM measurements are also easier to avoid than in the 
MSSM. One reason is that the complex phase in the effective \x parameter is 
dynamic. Hence it is possible that the phase is relatively small in the broken 
phase even though it varied strongly during the phase transition. Also, due 
to the stronger phase transition, electroweak baryogenesis is more efficient 
and the observed baryon asymmetry can be reproduced with smaller complex 
phases in the chargino sector. 

4.5.4 Summary 

In a probabilistic study, the collider and mass constraints provide quite 
strong bounds on the parameters of the scalar sector. However, once these 
constraints are passed, a large portion of the remaining parameter space 
leads to a strong first-order phase transition and viable baryogenesis in the 




(121) 
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nMSSM [104] . In this sense electroweak baryogenesis is a generic feature of 
the nMSSM. 

4.6 Other models 

For completeness we briefly mentioned in this section other models in which 
electroweak baryogenesis has been studied. This includes the Beyond MSSM 
scenario [!08[ll09j . the MSSM with an additional U(l)' gauge interaction pJIJJ 
11114 1112] . models with R-symmetric supersymmetry |113[ 1114] . the singlet 
Majoron model |115j and left-right symmetric models [116] . 
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5 Conclusions 

The main ingredients of electroweak baryogenesis are a strong first-order 
phase transition and new sources of CP violation. For this reason, elec- 
troweak baryogenesis is ruled out in the SM and heavily constrained in the 
MSSM. Nevertheless, in models with a more general scalar sector a strong 
first-order phase transition and electroweak baryogenesis are quite common 
features. 

From the perspective of electroweak baryogenesis, these models have the 
added benefit that the determination of the baryon asymmetry is much more 
robust than in the MSSM. In most of these models, the dominant source of 
CP violation arises from a semi-classical force that is sensitive to the spin of 
a single particle species. In contrast, in the MSSM the CP violation oper- 
ative during the phase transition arises from flavor mixing in the chargino, 
neutralino or stop sectors. This complicates the analysis through issues that 
are specific to systems with several flavors as flavor oscillations and resonant 
enhancements. 

Ultimately, whether electroweak baryogenesis is a realistic scenario hinges 
on the question if and how the hierarchy problem is solved by new physics 
at the electroweak scale. The LHC discovery of a Higgs-like particle of mass 
m = 125 GeV indicates that the MSSM can only solve the hierarchy problem 
at the cost of introducing a little hierarchy problem. This makes models with 
extended scalar sectors very attractive and in turn electroweak baryogenesis 
a promising mechanism for explaining the observed baryon asymmetry of the 
Universe. 
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A The weak sphaleron rate 

One essential ingredient of electroweak baryogenesis is the weak sphaleron 
rate 0. It couples to the left-handed fermions and anti-fermions of the Stan- 
dard Model and equally violates lepton and baryon number. In the presence 
of a (eventually local) CP-asymmetry in the left-handed particle densities, 
the sphaleron is biased towards a net baryon number. At the same time any 
pre-existing baryon number diffuses as long as baryon minus lepton number 
is conserved, B = L. The baryon asymmetry obeys the equation [29] 



where T ws is the weak sphaleron diffusion rate and denotes the chemical 
potential of the left-handed fermions. The final baryon asymmetry is then 
given by integration 



with g* ~ 106.75 the effective number of degrees of freedom at electroweak 
temperatures and we defined v = 45T WS /4:V W T 3 . The chemical potential \ll 
falls off at least as e~ DqZ in the symmetric phase where D q is the quark dif- 
fusion constant. Thus for large wall velocities v w , the exponent — vz is irrel- 
evant and the dependence on the wall velocity is inherited from the chemical 
potential \xl that is in leading order linear in v w . Hence for v <C D q and 
v w <C 1 the final baryon asymmetry depends only weakly on the wall ve- 
locity. If v w approaches the speed of sound, c s = l/\/3, diffusion should 
become inefficient (which however is not correctly reproduced in the analysis 
of sec. [3] that assumes small wall velocities). In the limit of very small wall 
velocities the exponent becomes important and leads to further suppression. 
This indicates that the wall is so slow that the sphaleron is saturated. In this 
regime backreactions on the left-handed chemical potential /il should not be 
neglected. 

On one hand, the sphaleron rate has to be large during the phase tran- 
sition in the symmetric phase in front of the wall. The CP violation in the 
reflection of particles leads to a net CP- violating particle density in front of 

8 An early review on the sphaleron rate in the context of electroweak baryogenesis is 
given in |117j . 
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the wall. If this particle density carries (positive) lepton or (negative) baryon 
number, the sphaleron process is biased towards a positive net baryon num- 
ber. In order to produce a baryon asymmetry of the observed magnitude 
7] ~ 10 -10 the sphaleron process should be considerably larger than 77 in 
electroweak units. 

This sphaleron rate in the symmetric phase was controversially discussed 
for some time in the literature 0. The main problem is that the sphaleron 
rate is non-perturbative due to the large occupation number of soft modes 
but also sensitive to the dynamics of the hard modes in the plasma [119] . The 
system is successfully described by Bodekers effective theory [120[ I121[ 1122] 
that can be easily simulated on a lattice. In conclusion, the weak sphaleron 
rate in the symmetric phase reads 

T ws = k f^-)a 5 w T\ (124) 

where m 2 D = ^-g^T 2 is the Debye mass of the weak gauge fields and g w is 
the gauge coupling of the weak interactions. Numerically the coefficient k is 
given by k — 40. Including the dynamics of the Higgs field slightly reduces 
this number and one finds [123J 

T ws ~ 1.0 x KT 6 T 4 . (125) 

This is in principle sufficiently fast for electroweak baryogenesis. 

On the other hand, the sphaleron rate in the broken phase should be 
smaller than in the symmetric one. For equal sphaleron rate no net baryon 
number would be generated during the phase transition, since the plasma 
in the bubble carries the opposite lepton and baryon number densities com- 
pared to the plasma in front of the wall. In fact, the sphaleron rate in the 
broken phase must be many orders of magnitude smaller than the rate in 
the symmetric phase. After the phase transition, the plasma components in- 
side and outside the Higgs bubbles mingle again. Even though a net baryon 
(and equal lepton) number was generated during the phase transition, the 
real equilibrium state of the system is still B = L = 0. If the sphaleron 
process is still active after the phase transition, the system has a time of 
order of the Hubble scale to attain this equilibrium. Hence, in order for elec- 
troweak baryogenesis to work, the sphaleron rate must be slow compared to 
the Hubble expansion. 



A nice summary of the status quo can be found in the talk [TT8] . 
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The sphaleron rate in the broken phase is accessible to semi-classical 
analysis [1241 1125] and is exponentially suppressed by the sphaleron energy 

T ws ~T 4 e-^ /T . (126) 

The sphaleron energy is proportion to [124] 

E, r ^~, (127) 
9w 

and numerically one finds H ~ 2.8. If one requires that the sphaleron rate 
is slow compared to Hubble expansion, T ws <C HT 3 , this leads to |125[ I126[ 
IT2T1IT28] 

0c>l.lT c , (128) 

This is the so-called sphaleron washout criteriorj^l. 

Also the sphaleron rate in the broken phase has been confirmed non- 
perturbatively on the lattice |131] . Recently, the first lattice calculations con- 
necting the symmetric phase with the broken phase have been presented [132J, 
confirming the picture developed in the two different phases in a unifying 
framework. 



B Semi-classical approach to phase transitions 



The formalism to describe semi-classical tunneling was pioneered in con- 
densed matter systems by Langer [133], in quantum field theory by Coleman 
|134[ 1135] and at finite temperature by Linde |136j . A review of the topic 
can be found in |137] . 

In a tunneling problem the effective potential has at least two local min- 
ima that constitute the different phases the physical system can reside in. 
In the following we call these two phases the symmetric (before the phase 
transition) and broken (after the phase transition) phases, motivated by the 
electroweak phase transition, see Fig[12j In the semi-classical WKB approx- 
imation, the tunnel probability per volume and time is suppressed by the 
Euclidean action of the so-called tunneling bounce 

P-Ae"^, (129) 
10 For a more detailed discussion of this argument see also |129[ 1130] . 
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Figure 12: Example for a potential with a metastable minimum. The phase 
transition proceeds from the symmetric phase <f) s to the broken phase fa. 

derived from the effective action expanded in gradients 

S ~ J d A x ^06^0 + V{<j>) , (130) 

where V(4>) denotes the corresponding effective potential that eventually de- 
pends on the temperature. 

The coefficient A in (I129P is for dimensional reasons of electroweak scale, 
A ~ T 4 . The phase transition happens when the probability to nucleate a 
bubble of broken phase is of order unity in a Hubble volume and time leading 
to the condition 

S -log A -140. (131) 

The bounce is at zero temperature a 0(4)-symmetric solution to the Eu- 
clidean equations of motion while at finite temperature it is 0(3) symmetric 
and periodic in imaginary time. The equations of motion then read 

d 2 <j> | , d ^ d4 _ dV 
dr 2 t dr dtp ' 

with d — 4 (d — 3) for tunneling at zero (finite) temperature. The bound- 
ary conditions are such that <j) initially rests close to the broken phase and 
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asymptotically approaches the symmetric phase at late 'time' 

0(0) 0'(O) = O, 0(oo)=0 s . (133) 

In the limit of weak phase transitions, the thin-wall approximation ap- 
plies |134j . In this case the field <p rests for a rather long time tr close to 
the broken phase and then quickly changes to the symmetric phase. In this 
case, the tunnel action can be reexpressed in terms of the wall tension 



a 



d<P^2V(<P) , (134) 



and the potential difference AV = V(<pb) — V{<j>s) as 

s = W < d = 3 >' (135) 

Otherwise, for one field and quite arbitrary conditions, the tunneling action 
can easily be obtained numerically using the shooting-algorithms [134] . For 
several scalar fields, more involved methods have to be used [1381 1139j . 

Recently, the gauge-independence of above approach was questioned [1291 
1140^ H41j but an explicit calculation in a Abelian toy model shows that the 
dependence on the gauge choice is actually quite small [142] . This is also 
supported by the fact that the semi-classical approach agrees reasonably well 
with non-perturbative methods on the lattice [143] . The main corrections to 
the procedure above seem to come from higher order contributions to the 
kinetic term and the effective potential in the effective action f)130p . 



C Wall velocity and wall thickness 

Several parameters of the phase transition enter the produced baryon asym- 
metry quantitatively. Namely the critical vev <p c , the critical temperature T c , 
the wall thickness £ w and the wall velocity v w . 

The most important one is hereby the ratio <p c /T c that determines the 
sphaleron washout and also the reflection of the particles by the Higgs wall 
that in turn leads to the CP violation in the particle densities. Fortunately, 
in most models with viable electroweak baryogenesis this quantities can be 
rather easily obtained using the semi-classical methods of appendix [Bj 
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Another important input is the wall thickness £ w . The gradient expansion 
can only be applied for thick walls, £ w T ^> 1, and the final baryon asymmetry 
is in the one-flavor case roughly inversely proportional to the wall thickness. 
For not too large wall velocities, the shape of the Higgs bubble profile does 
not change much during the expansion [144] . The wall thickness can then 
be determined from the wall thickness of the nucleated bubbles in the semi- 
classical tunneling analysis. 

Finally, the wall velocity v w enters the analysis. Viable baryogenesis re- 
quires that the wall velocity is small enough to allow for particle diffusion in 
front on the wall v w < 1/ \/3. For wall velocities smaller than that the pro- 
duced baryon asymmetry is rather insensitive to the wall velocity as already 
discussed in appendix |A] This results from the fact that the CP violation 
accumulated in front of the wall is proportional to the wall velocity. At the 
same time, the phase transition proceeds longer and the sphaleron process 
can act longer on the CP-violating particle densities and convert them into a 
baryon asymmetry. In this regime the final baryon asymmetry depends only 
weakly the wall velocity. However, for very slow walls, the sphaleron process 
becomes saturated and the final asymmetry scales linearly with the wall ve- 
locity. Due to the smallness of the sphaleron rate, this typically happens for 
velocities of order v w < 10 -3 . 

So the pivotal question is if the wall velocity is in the regime 10~ 3 <C v w < 
l/y/3 where above approximations are reasonable and the final asymmetry 
is insensitive to the wall velocity. To answer this question in a specific model 
requires to perform an out-of-equilibrium analysis that so far was only per- 
formed in the Standard Model |145[ |4"8] and the MSSM [146] . In both cases, 
the wall velocity turned out to be in the desired ballpark. For other models, 
the wall velocity is still unknown. A simple way of estimating the wall veloc- 
ity is to model friction in a phenomenological approach and to extrapolate 
the results from the SM and the MSSM [HO EHJ Ell [T29lfT50] . 

D Electroweak phase transition in the SM 

In this section we review the perturbative analysis of the electroweak phase 
transition in the Standard Model. We follow the work [151] but present a 
simplified analysis. 
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At tree level, the effective potential of the Higgs field reads 

V = ±(<p 2 -v 2 ) 2 , (136) 
and at one loop order the thermal corrections to the free energy are 

AV 1 = Tt^2 / dxx2 log ( X ± ex P(--\/ x2 + m f^ , ( 137 ) 

where ± stand for fermions/bosons respectively, T denotes the temperature, 
j3 the inverse temperature and mj the different particle masses. As long as 
the masses do not exceed the temperature, this can be expanded as 

—m 2 T 2 + 0(m 4 ) , 
48 v ' ' 

—m 2 T 2 - — m 3 T + 0(m*). (138) 
24 12tt V ' K ' 

Of special importance are hereby the cubic terms contributed by the bosons. 
If the mass of a bosonic field is only generated by the coupling to the Higgs 
vev (as is the case for the weak gauge bosons in the SM), this gives in turn 
rise to a term of the form (p 3 T in the effective potential. This term is essential 
to generate a potential barrier between the symmetric and the broken phase. 
Consider a potential of the form 

V = /j, 2 (T)(f) 2 -ET(f) 3 + j(j) A . (139) 

At some temperature T c , this polynomial potential has two degenerate min- 
ima at <p = and <fi = <ft c > and is of the form 

V=^4> 2 (<P-<p c ) 2 . (140) 

Comparison with (11391) then shows that 

/i 2 (T c ) = iA^, ET c = h<j> c . (141) 



AV} 

* bosons 



This immediately implies 

C /T C = 2E/X (142) 
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and larger Higgs masses lead to weaker phase transitions. 

In the Standard Model, the cubic coefficient arises only from the (trans- 
verse) electroweak gauge bosons [152[ 1153] 1154] . E ~ 1CT 2 . Accordingly, a 
phase transition strong enough for electroweak baryogenesis is only possible 
for Higgs masses below 40 GeV [155] in light of the constraint (I128p . Besides, 
for Higgs masses rrih > 70 GeV the perturbative analysis breaks down and a 
cross-over replaces the phase transition. 

In the MSSM, additional contributions to the cubic term come from the 
right-handed stops in case their mass is below the top mass. This can make 
the phase transition strong enough for baryogenesis for some parts of the 
parameter space, even when a Higgs of mass ~ 125 GeV is assumed (see 
sec. H3J. 

In general, if the potential barrier arises from a thermal cubic contribu- 
tion, the relation ( 1142 j) in combination with a Higgs mass of rrih — 125 GeV 
implies that a strong first-order phase transition requires at least E > 0.1. 
So a moderate number of light bosons that couple strongly to the Higgs is in 
this case essential. Yet, in many models the strength of the phase transition 
does not rely on the thermal cubic contributions. The prime example for this 
are models with an extended Higgs sector. When several scalar fields acquire 
a vev at electroweak scales, potential barriers can arise even in the tree level 
scalar potential (see sec. I4.2p . 
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